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ABSTRACT 


AN EXPERIMENTAL AND ANALYTICAL METHOD FOR 
APPROXIMATE DETERMINATION OF THE TILT 
ROTOR RESEARCH AIRCRAFT 
ROTOR/WING DOWNLOAD 

David E. Jordan 
May 1986 

An investigation was performed to analyze the XV-15 
Tilt Rotor Research Aircraft download phenomenon. This 
phenomenon is a direct result of the two rotor wakes 
impinging on the wing upper surface when the aircraft is in 
the hover configuration. This study represents one of many 
devoted to the analysis of this problem. 

For this study the analysis proceeded along two lines. 
First was a method whereby results from actual hover tests 
of the XV-15 aircraft were combined with drag coefficient 
results from wind tunnel tests of a wing that was represen- 
tative of the aircraft wing. Second, an analytical method 
was used that modeled the airflow caused by the two rotors. 
Formulas were developed in such a way that a computer 
program could be used to calculate the axial velocities in 
the vicinity of the aircraft's wing. The axial velocities 
were then used in conjunction with the aforementioned wind 
tunnel drag coefficient results to produce download values. 
An attempt was made to validate the analytical results by 
modeling a model rotor system for which direct download 
values were determinable. 



Table of Contents 


Page 

List of Tables v 

List of Figures V1 

List of Symbols viii 

Chapter 

1. Introduction 1 

Description 1 

History 4 

2. Experimental Determination of the Download . 9 

Hover Test 9 

Wind Tunnel Tests 20 

Experimental Download Analysis 22 

3. Analytic Determination of the Download ... 28 

Vortex Ring Induced Velocity 29 

Wake Contraction 33 

Time Spacing of Vortex Rings 35 

Analytic Download Analysis 36 

Computer Program Calibration 39 

Download Computer Program 41 

Velocity Profiles 44 

4. Model Test 49 

5. Results and Discussion 51 

Experimental Download Analysis 51 


i i i 



Chapter Page 

5. (continued) 

Model Test Analysis 55 

Analytic Download Analysis 56 

6. Conclusion 63 

References 66 

Appendixes 

A. Transformation from Equation 6 to 

Equation 7 68 

B. Empirical Adjustment to Power Polar to 

Reflect Nonlinear Aerodynamics Region ... 70 

C. Development of the Biot-Savart Equation 

to Determine the Axial Velocity 

Component Due to a Vortex Ring 73 

D. Calibration Version of XV-15 Computer 

Program 82 

E. Download Version of XV-15 Computer 

Program 88 

F. Calibration Version of Model Rotor 

Computer Program 100 

G. Download Version of Model Rotor 

Computer Program 106 


iv 



Tables 


Table Page 

1. Out of Ground Effect Full Span Flap 

Hover Data 12 

3/2 

2 . Deviation of C From Linear C Vs . C 

Relationship w p . . . 7 . . . . 20 

3. Rotor Lifted Weights and Flap Deflection 

Angles at a Constant Power Setting .... 21 

4. Drag Coefficients and Flap Deflection 

Angles 23 

5. Experimental Download 24 

6. Analytic Download 58 


V 



Figures 


Figure Page 

1. The XV-15 Tilt Rotor Research Aircraft ... 2 

2. XV-15 Flap/Flaperon Schedule 10 

3. XV-15 Engine Exhaust Jet Thrust 13 

4. XV-15 Power Coefficient /Weight 

Coefficient Polar 18 

5. Weight Coefficient Deviation From Linear 

Aerodynamics as a Function of Power 
Coefficient 19 

6. XV-15 Experimental Power Coefficient/ 

Weight Coefficient Polar 19 

7. Effect of Flap Angle Deflection on 

Rotor Lifted Weight 21 

8. Effect of Flap Angle Deflection on Two- 

Dimensional Drag Coefficient 23 


9. Effect of Flap Angle on Drag Coefficient 
Change, Referenced to the Drag 

Coefficient for 0° Flap Deflection .... 25 


10. Change in Rotor Lifted Weight, Referenced 

to Rotor Lifted Weight for 0° Flap 

Deflection, as a Function of the 

Change in Drag Coefficient 26 

11. Experimental XV-15 Download As a Function 

of Flap Deflection Angle 27 

12. The Geometry Used for Deriving the Axial 

Velocity Component Due to a 

Vortex Ring 30 

13. Vortex Ring Downstream Position as a 

Function of Time 38 

14 . In Ground Effect Vortex Ring 

Representation . 43 

15. Computer -Generated Velocity Profile, 

C t = .0094 46 


- VI 



Figure 


Page 


16. Computer-Generated Velocity Profile, 

C t = .01 47 

17. Computer-Generated Velocity Profile, 

C t = .0105 48 

18. Model Rotor Test Setup 49 

19. Effect of Hover Height on Analytic 

XV-15 Download 58 


vi i 



Symbols 



2 

rotor disk area (ft. ) 
drag coefficient 
power coefficient 
profile power coefficient 


C^: thrust coefficient 

C w : weight coefficient 

c: wing chord (ft.) 

DL: aircraft rotor /wing download (lb.) 

DL/T: aircraft download-to-thrust ratio 

4^: flap deflection angle (deg.) 

Fjj engine exhaust jet thrust (lb.) 

GW: aircraft gross weight (lb.) 

2 

r : circulation strength (ft. /sec.) 

H/D: rotor height above ground to rotor diameter ratio 

IGE: in ground effect hover 

K: vortex ring time spacing factor 

OGE: out of ground effect hover 

R: rotor radius (ft.) 

r': vortex ring radius (ft.) 

3 

p : ambient air density (slugs/ft. ) 

STOL: short take-off and landing 

T: rotor thrust (lb.) 

VTOL: vertical take-off and landing 
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Symbols (continued) 

average induced velocity at the rotor disk 
( ft ./sec . ) 

rotor tip speed (ft. /sec.) 
aircraft rotor lifted weight (lb.) 



CHAPTER 1 


Introduction 

Description 

The XV-15 (Figure 1) is a "Tilt-Rotor" aircraft. The 
aircraft has a wingspan of 32 feet 2 inches, a length of 46 
feet 3 inches, and a design gross weight of 13,000 pounds. 
The aircraft's propulsion system consists of two Lycoming 
LTC1K-4K turboprop engines and two 25-foot diameter, three- 
bladed prop-rotors. These engines are variants of Lycoming 
T53-L-13B engines and are mounted one on each wing tip. 

The rotor axes are capable of rotation from a vertical 
position (for helicopter and hover flight) to a horizontal 
position (for airplane flight). Control of the aircraft in 
the helicopter configuration is provided by rotor-generated 
forces and moments, while in the airplane configuration 
control is provided primarily by conventional control sur- 
faces. Operation at engine tilt angles between vertical 
and horizontal provides short take-off and landing (STOL) 
capabilities, as opposed to vertical take-off and landing 
(VTOL) capabilities. Cross shafting of the two engines 
provides rotational speed synchronization, permits power 
transfer for transient conditions, and prevents complete 
loss of power to either rotor due to a single engine 
failure . 
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The XV-15 was designed and built as a "proof of 
concept" aircraft and is providing a data base for the 
larger scale V-22 Osprey that is currently in the prelimi- 
nary design phase. Two XV-15 aircraft were built by Bell 
Helicopter Textron Inc. under contract by NASA. One, air- 
craft 702, was retained by Bell Helicopter and is operated 
in Texas; the other, aircraft 703, is operated by NASA at 
Ames Research Center. 

The XV-15 tilt rotor research aircraft is a unique 
aircraft with unique capabilities. Along with the capabil- 
ities, however, come unique problems. One problem occurs 
when the aircraft is operating in a hover. In this con- 
figuration the engines are in a vertical position with the 
rotor disks parallel to the wing surface. The wing is 
approximately 5 feet below the rotor disk and experiences a 
"downward" directed force due to the direct impingement of 
the rotor wake on the wing. The approximate determination 
of the magnitude of this downward directed force, hereafter 
referred to as download, was the goal of this research 
project. The download and its elimination, or at least 
reduction, is a significant problem for today's XV-15 proof 
of concept prototype and tomorrow's V-22 Osprey production 
aircraft. 

It is difficult to experimentally determine the air- 
craft download due to the number of variables involved: 
rotor thrust, engine exhaust jet thrust, aircraft gross 
weight, download, and upload in the case of a ground effect 



reflected wake. The problem with the experimental deter- 
mination of the download arises from the fact that there is 
no readily apparent way to separate the download out of a 
gross measurement of the vertical force system. While the 
aircraft gross weight can be exactly determined, the 
exhaust jet thrust can only be approximated, being a 
function of inlet air velocity, and the installed rotor 
thrust defies determination due to its direct interrelation 
with the download phenomenon. These facts have led to many 
attempts to determine download through model test and theo- 
retical techniques. 


History 

The first record of a download type phenomenon was 
encountered in the late 1950s with the XV-3 (Deckert and 
Ferry, 1960). The XV-3 was a S/VTOL aircraft much like the 
XV-15 except that it was powered by one ra’dial type piston 
engine mounted in the aircraft fuselage. Power was shafted 
to a transmission at each wing tip which drove two bladed 
prop -rotors . 

During the flight test evaluations of the XV-3 it was 
discovered that a large increase in power required was 
encountered suddenly when hover flight was approached. 

This fact was attributed to an aspect ratio effect. The 
reasoning was that theoretically, above transitional 
speeds, the prop-rotors were actually influencing a 
diameter of air extending from the outboard edge of one 
prop-rotor to the outboard edge of the other. With this 



large mass of air being influenced, the induced power 
required was necessarily low. However, below transitional 
speed (i.e., in hovering flight), each prop-rotor was 
thought to influence a diameter of air approximately equal 
to its own diameter. Thus, with the much smaller mass of 
air being influenced, the induced power requirement had to 
be increased to produce the same thrust. While this expla- 
nation seems feasible, the problem could have been the 
first encounter with tilt-rotor download. 

In 1969 a test was conducted at Wright-Patterson Air 
Force Base that placed a mock-up under a full-scale CH-47 
Chinook rotor (Gillmore, 1971). It was speculated that on a 
wing without leading or trailing edge devices the download 
force would approach 12% of the aircraft gross thrust. The 
test wing, however, was run utilizing trailing edge flaps 
deflected 70° in conjunction with a leading edge device that 
opened like an umbrella, top and bottom. Using this con- 
figuration it was demonstrated that the download could be 
reduced to 4% of gross thrust. While these results seem 
perfectly accurate and acceptable, the extension of these 
values to reflect the XV-15 would be questionable. The 
test used only one rotor and no image plane or end plate to 
simulate the longitudinal symmetry of the XV-15. Thus, the 
test did not adequately model the flow phenomenon that 
exists on the XV-15. 

Other tests were conducted by the Bell Helicopter 
Company beginning in 1968 (Marr, Sambell , and Neal, 1973) 



that utilized a 1/5 scale powered aeroelastic model of the 
XV-15. The model was tested at rotor height above ground 
to rotor diameter (H/D) ratios from 0.5 through 1.0. In an 
attempt to isolate the download force, the model was run 
both with and without wing panels installed. The test 
revealed that at an H/D value of 1.0 the wing panel down- 
load was 11.4% of the model thrust. This value was higher 
than the 7% that was predicted for the full-scale aircraft. 
It was thought that the difference between the predicted 
and actual values could be attributed to a Reynold's number 
effect. A Reynold's number correction had the effect of 
reducing the predicted download value to 5-7% of the full- 
scale thrust. The results of these tests are also suspect 
because of two factors. First, the Reynold's number cor- 
rection was based on wind tunnel data of "cross sectional 
shapes with rounded edges and flat sides" and not the model 
aircraft wing. Second, the download predicted for the 
full-scale tilt rotor aircraft was based on in ground 
effect (IGE) data and therefore does not adequately reflect 
the general out of ground effect (OGE) solution. 

The latest tests that addressed the download issue 
were performed at the Outdoor Aerodynamic Research Facility 
(OARF) at NASA-Ames Research Center in 1984 (McVeigh, 

1985). The download aspect of the tests was specifically 
directed toward the V-22 Osprey program but it was felt 
that the results were applicable to this study. The tests 
consisted of a 25 -foot diameter three-bladed prop-rotor 



mounted with its axis of rotation parallel to the ground. 
This configuration precluded any influence that the ground 
would have had on the prop-rotor performance. The effect 
of the opposing prpp-rotor was simulated by using a large 
image plane. A wing, complete with a 31% of wing chord 
full span flap and follower, was placed downstream of the 
rotor disk between the rotor axis and the image plane. The 
wing was instrumented to facilitate the measurement of 
forces and moments . The results of the test show that for 
a flap deflection angle of 67° and a nominal design thrust 
coefficient of .016, the download-to-thrust ratio was 9.3%. 
The wing used for the tests had a chord of 10.5 feet and 
was placed approximately 8 feet downstream of the rotor 
disk with a 6° dihedral. This test represents the most 
recent and most accurate work done to date attempting to 
resolve the download issue. 

The above tests show that the download problem is 
significant and deserves attention. The fact that early 
estimates of the download were in the 3% to 5% range but 
have grown to 9% to 10% of gross thrust is significant in 
itself. Applied to a hovering XV-15 aircraft, the dif- 
ference in the early and late estimates of download account 
for a difference of almost 1,000 pounds. 

The download analysis, for this study, was performed 
utilizing two different procedures. The first method was 
an experimental method that combined results from actual 
aircraft hover tests of the XV-15 and results of wind 



tunnel tests of a model XV-15 wing. The second method was 
a computer simulation of the aircraft propulsion system, 
namely two 25-foot diameter prop-rotors, that relied on 
vortex ring theory. This method also made use of the same 
wind tunnel test results. 



CHAPTER 2 


Experimental Determination of the Download 

Hover Test 

Flight tests were utilized to attempt to experimen- 
tally determine the XV-15 rotor/wing download. A hover 
test was conducted in 1983 in which the aircraft was flown 
using five different deflection settings with careful 
attention to maintain the aircraft gross weight within a 
narrow range (approximately 13,500 to 13,700 pounds). Data 
from the hover test was utilized in conjunction with data 
acquired from wind tunnel tests of a representative two- 
dimensional airfoil to arrive at values of download-to- 
thrust ratio. 

The XV-15 aircraft wing is equipped with a flap/ 
flaperon control surface system. The system is rigged to 
provide a linear relationship between the flap and flaperon 
deflection, i.e., a flap deflection of 40° causes a flaperon 
deflection of 25°, a flap deflection of 20° causes a flaperon 
deflection of 13°, etc. (Figure 2). Throughout this linear 
deflection relationship, roll control is accomplished by 
the pilot's flaperon input overriding the associated flap/ 
flaperon angle relationship. The flaperon deflection, 
corresponding to the pilot input, is from the position 
predetermined by the flap setting. For the hover download 
investigation, a special modification was implemented to 
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the flap/flaperon system which facilitated ground adjust- 
ment of the flaps and flaperons to the same deflection 
angle, thereby producing the effect of a full span flap 
system (Figure 2). The deflection angles used during the 
test for the full span flap system were 0°, 20°, 40°, 50°, 
and 60°. 



Flap Angle £ t - 

(deg.) 

Figure 2 

XV-15 Flap/Flaperon Schedule 

In order to facilitate a direct correlation between 
test points throughout the hover test, efforts were made to 
keep the aircraft gross weight within a narrow range. 
Ballast was added to the aircraft during the flap readjust- 
ment ground time to compensate for the fuel burned during 
each portion of the test. Accordingly, the aircraft gross 
weight was maintained within a 200-pound range. Aircraft 
gross weight determination during the test sequences was 



performed after the tests using on-board fuel quantity 
instrumentation data. 

An aircraft wheel hover height of 50 feet above the 
ground was established for the tests using visual 
references. This height corresponds to an H/D ratio of 
2.5, insuring that the data taken would represent OGE 
conditions . 

Rotor speeds were held constant (+1 RPM.) throughout 
the test points through the use of the engine RPM governor. 
Three data points per flap setting were recorded and are 
presented in Table 1. Ambient test conditions were moni- 
tored and recorded through the use of a local weather 
station that was placed well outside of the rotor influ- 
ence. Light winds existed throughout the test period 
ranging from 2.5 knots at 320/045° to 4 knots at 340°. The 
aircraft heading was held constant at 320° during the tests. 

The system of vertical forces that act on the XV-15 
during hover operation is composed of four components: 
engine exhaust jet thrust, rotor thrust, aircraft gross 
weight, and the rotor/wing download. 

Engine jet thrust was determined with an acceptable 
degree of accuracy using Lycoming LTC1K-4K engine data 
(Figure 3). Engine exhaust jet thrust is dependent on 
engine horsepower and free stream velocity in the vicinity 
of the engine inlet. After the rotor tip speed and mast 
torque were determined from aircraft instrumentation, rotor 
horsepower could be calculated. Assuming that the 



Table 1 


Out of Ground Effect Full Span 
Flap Hover Data 


RECORD 

FLAP 

(DBG) 

GROSS WT 
(LBS) 

WIND 

(V k /DJRECT) 

o 

ROTOR 

SPEED 

(RPM) 

MAST TORQUE (WLB) 
LEFT RIGHT 

534 

0 

13550 

2.5/045° 

0.9917 

587 

116814 

118841 

535 

0 

13536 

2.5/045° 

.9917 

587 

114808 

116764 

536 

0 

13527 

2.5/045° 

.9917 

587 

115728 

118494 

537 

20 

13537 

3.5/340° 

.9917 

587 

112384 

114652 

538 

20 

13521 

3.5/340° 

.9917 

587 

113043 

114092 

539 

20 

13513 

3.5/340° 

.9917 

587 

114356 

115005 

54 0 

40 

13657 

3.0/345° 

.9955 

587 

111464 

113951 

541 

40 

13647 

3.0/345° 

.9955 

587 

111305 

114479 

54 2 

40 

13637 

3.0/345° 

.9955 

587 

110827 

115429 

54 3 

50 

13673 

3.0/340' 

.9917 

587 

109476 

112615 

544 

50 

13662 

3.0/340° 

.9917 

587 

108636 

111456 

545 

50 

13645 

3.0/340° 

.9917 

58 7 

109314 

113052 

547 

60 

13614 

4.0/340° 

.9936 

588 

107176 

118768 

548 

60 

13601 

4 .0/340° 

.9936 

588 

108460 

108746 

549 

60 

13585 

4.0/340' 

.9936 

587 

105907 

109352 


Parameters 


RECORD 

AVG. 

SHAFT 

POWER 

(HP) 

c p 

C P 

ENGINE 

JET 

THRUST 

(LBS) 

WEIGHT 
(GW-2F j ) 

C w 

(NON 

LINEAR) 

C w 3/2 

C v 

(LINEAR ) 

534 

1098 

.001150 


85.0 

13380 

.00979 



535 

1079 

.001130 

.001141 

83.5 

13369 

.00978 

.000967 

.01007 

536 

1091 

.001143 


84.0 

13359 

.00977 



537 

1057 

.001108 


83.0 

13371 

.00978 



538 

1058 

.001108 

.001112 

83.0 

13355 

.00977 

.000966 

.010043 

539 

1068 

.001119 


83.0 

13347 

.00977 



540 

1050 

.001095 


82.5 

13492 

.00983 



541 

1052 

.001097 

.001097 

83.0 

13481 

.00982 

.000973 

.01006 

542 

1054 

.001099 


83.0 

13471 

.00981 



543 

1034 

.001083 


82.0 

13509 

.00988 



544 

1025 

.001074 

.001080 

81.7 

13499 

.00989 

.000982 

.010107 

545 

1036 

.001084 


82.0 

13481 

.00986 



547 

1008 

.001048 


81.0 

13454 

.00979 



548 

1013 

.001054 

.001050 

81.5 

13438 

.00978 

.000969 

.009987 

549 

1003 

.001048 


80.2 

13425 

.00980 





200 - 


V =0 knots 


V' 0 = Engine Inlet Free Stream 
Velocity 

Data Source: Lycoming LTC 1K-4K 
Model Specification Number 104.47 
Figure 1.6.1 


Figure 3 

XV-15 Engine Exhaust Jet Thrust 


transmission linking the engine and rotor was 93% effi- 
cient, calculation of the engine horsepower was then 
possible. Computation of the free stream velocity in the 
vicinity of the engine inlet was accomplished through the 
use of simple momentum theory (i.e., v^ = V. C^/2 ) . 

In order to determine the aircraft thrust coefficient, the 
thrust per rotor was assumed to be 110% of half the air- 
craft gross weight and the rotor tip speed was determined 
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from aircraft instrumentation. These two values facili- 
tated calculation of the average rotor-induced velocity. 
This value of average rotor-induced velocity was assumed to 
approximately represent the free stream velocity in the 
vicinity of the engine inlet. With the induced velocity 
and engine horsepower determined, a first approximation of 
engine exhaust jet thrust was possible. Further iterations 
to refine the engine exhaust jet thrust values were deemed 
unnecessary due to the apparent insensitivity of the jet 
thrust to the level of rotor thrust and also the relatively 
small value that the jet thrust represented compared to the 
other forces in the system. Engine jet thrust during the 
hover tests was determined to be between 80 and 85 pounds 
per engine. Doubling this value accounted for the total 
engine exhaust jet thrust. This total engine exhaust jet 
thrust was then subtracted from the aircraft gross weight 
and the value that resulted represented the "rotor lifted" 
aircraft weight (i.e., W = GW - 2Fj). 

The rotor lifted aircraft weight was used to define a 
"weight coefficient." This weight coefficient was in the 
form of the standard helicopter thrust coefficient, with 
half of the total rotor lifted weight substituted for the 
rotor thrust. The form of the weight coefficient was 

c ■ 

" 2 ’V? 


( 1 ) 
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In order to be able to apply the OGE hover test data 

3/2 

at other levels of power coefficient (C^) , the Cp = mC t 

+ C power polar algorithm (Hayden, 1976), applicable only 
p o 

in the linear aerodynamic range, was utilized. This 
algorithm expresses the total power coefficient as being 
composed of the induced power, which is proportional to the 
thrust coefficient to the three halves power, and the 
profile power, which is that power required to turn the 
rotor in a zero thrust condition. Although this is an 
approximation, for hovering flight the combination of 
induced and profile power accounts for 90% to 100% of the 
total power. 

Solving the algorithm for thrust coefficient (C. ) 


gi ves : 


S - V 


Calculating the slope (m) between a general point 
„ 3/2 . . 

^ K rv P 4 t-\ ♦- f P \ i «»a r- 


(C P ' C t 
P 1 ^1 


) and the C intercept (0, C ) gives: 

P p 3 

o 


C - C. 


Substituting this expression into that for gives: 


C - C 


(C t ) 
l l 


where (C , C. ) is a general OGE hover point corresponding 
P 1 *1 

to a given flap setting and C is the C axis intercept of 

P Q P 
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3/2 

the = mC t + C curve, which is based on prior XV-15 
hover data (Hoerner, 1965). Referring back to the 
relationship that the rotor lifted weight is equal to the 
aircraft gross weight less the engine exhaust jet thrust 
plus the rotor/wing download gives 

T = W + 01 ( 5a ) 

or 

DL 


T = M(l + -{H 


(5b) 


Dividing through by the density, rotor disk area and rotor 

2 

tip speed squared ( pApV t ) gives: 


V ‘J. ♦ tH 

Using the manipulation outlined in Appendix A and 
substituting this final result into the previous 
relationship for thrust coefficient (C^) gives: 

0L 


( 6 ) 


C = C (l - ^ |~* 

t t J 


(7) 


Due to the fact that it has been shown (Makofski and 
Menkick, 1956) that the rotor/wing download to rotor thrust 
(DL/T) ratio is independent of the rotor disk loading but 
dependent on body configuration in the rotor wake alone, 
this ratio for the XV-15 is assumed to be dependent on flap 
configuration only. Thus (1 + DL/T) was assumed to have a 
constant value for each flap angle setting and the expres- 
sion for C fc derived from the power polar expression became: 


C P - C P 0 

C_ - C 


2/3 


(C w > 


( 8 ) 



This relationship, depicted in Figure 4, states that the 

rotor lifted aircraft weight (W) is, for a specific flap 

setting, only a function of applied power, Cp. Therefore, 

for each flap setting and its associated C and C 

P 1 W 1 

values, a unique value will be defined. From this it 

follows that for each unique C w there will be a unique W, 

or rotor lifted aircraft weight. This relationship is 

particularly useful due to the fact that weight coefficient 

(C ) values were easily obtainable from the hover test data 

(Table 1) while true thrust coefficient values could not be 

initially determined. The linear relationship of the 
3/2 

C = mC fc + C expression is applicable only for values 
p p o 

of C below .0008, as can be seen in Figure 5. The 
P 

3/2 

deviation of C from the linear C vs. C relationship 

w p w 

is presented in Figure 6. At values of Cp above .0008, 
which is the case with the OGE hover data taken for this 
study, an empirical adjustment was utilized to reflect the 
departure from linear aerodynamics (Table 2). This adjust- 
ment consisted of finding a Cp dependent 4C w value from 
Figure 6, utilizing the study values of Cp. This a C w was 

then added to the nonlinear C values to arrive at the 

w 

3/2 

associated C on the C vs . C linear projection. A new 

w p w J 

value of C w , at the desired Cp, is then obtained by inter- 
polation (using the linear equation. Equation (8)). A 

Becond adjustment was then made to project the interpolated 

3/2 

C back to the nonlinear C vs . C curve. The second 

w p w 

AC value was obtained from Figure 6 at the desired C 
w p 














Table 2 


Deviation of G From Linear C Vs. 

Relationship ^ 



C P 

c 3/2 

w 

UN 

r 3/2 
Sr 

ACT 

c w 

UN 

C w 

ACT 

C w 

.0008 

.0007 

.0007 

.00788 

.00788 

0 

.0009 

.000808 

.000796 

.00867 

.00859 

.00008 

.0010 

.000908 

.000888 

.00938 

.00924 

.00014 

.0011 

.001012 

.000972 

.01008 

.00981 

.00027 

.0012 

.001114 

.001055 

.01075 

.01036 

.00038 

.0013 

.001215 

.001134 

.01139 

.01087 

.00051 

.0014 

.001316 

.001212 

.01201 

.01137 

.00064 

.0015 

.001418 

.001283 

.01262 

.01181 

.00081 

( .001) . 

Utilizing 

this method. 

outlined 

in Appendix 

B, 


computed values of C w and W were determined and are 
presented, with the associated flap angle test points, in 
Table 3. The incremental weight (aW) column in Table 3 
reflects the additional weight, referenced to the 0° flap 
angle weight, that could be hovered OGE as a result of 
increased flap angle deflection. The effect of flap angle 
on rotor lifted weight (aW) is presented in Figure 7. 


Wind Tunnel Tests 

Wind tunnel tests of two-dimensional wing sections 
were performed in the Army 7-foot by 10-foot wind tunnel at 
NASA Ames Research Center during April and May of 1983. 

The wind sections utilized plain flap configurations and 
were tested at angles of attack near -90°. Flaps of various 
chord sections were tested. 








The variation of drag coefficient with flap angle for 
the various flap chords tested is presented in Figure 8. 

The decrease in drag coefficient with increasing flap 
deflection angle is interrupted by a sharp rise, reflecting 
the onset of separated flow on the flap upper surface. The 
wing section with the .25c flap, which is the same size 
flap as that on the XV-15, demonstrated flow separation 
characteristics about 10° earlier than observed on the air- 
craft. The cause of this could be traceable to the fact 
that the XV-15 flap is slightly thicker than that which 
would be dictated by the basic airfoil contour. The result 
of the greater flap thickness is a larger flap leading edge 
radius than that which was built into the wind tunnel flap. 
Since the drag coefficient characteristics for each of the 
flap chords tested were similar in the attached flow 
region, it was assumed that the XV-15 wing drag coeffi- 
cients could be represented by extrapolation of the .25c 
flap data to 60° flap deflection, based on the performance 
characteristics of the ,30c and .35c flaps. The extrapola- 
tion is shown as a broken line in Figure 8 and is included 
in Table 4. 


Experimental Download Analysis 
The AC^ values listed in Table 4 represent drag 
coefficient decreases, referenced to the 0° flap deflection 
drag coefficient, for each increase in flap deflection. 

The relationship of AC^ to flap deflection angles is 



Flap Angle St 
(Jog.) 


Figure 8 

Effect of Flap Angle Deflection on Two 
Dimensional Drag Coefficient 

Table 4 

Drag Coefficients and Flap 
Deflection Angles 







depicted in Figure 9. Realizing that the decrease in drag 
coefficient ( AC^) was responsible for the increase in 
rotor lifted weight ( AW) for a given power setting (C^) , a 
relationship between AC^ and AW was established and is 
shown in Figure 10. Using a linear regression technique, 
it was possible to establish a ( AW/ a ) ratio, which was 
valid and constant for a specific power setting. Multipli- 
cation of this ratio by the two-dimensional wing drag 
coefficient for each associated flap deflection angle made 
it possible to define a form of the rotor/wing download as: 


DL 





(9) 


Results of the experimental download analysis for the 
selected power setting (Cp = .001) were computed as 
percentages of the total aircraft thrust, which was a 
function of flap deflection angle. The computed download 
values are listed in Table 5 and depicted in Figure 11. 


Table 5 

Experimental Download 


4 f W DOWNLOAD TOTAL THRUST DOWNLOAD 

(GW-2F. ) (AW/AC. )C. (W+DL) (% THRUST) 

3 a d Sf 


0 

12402 

1885 

14 287 

13.2 

20 

12599 

1706 

14305 

11.9 

40 

12753 

1472 

14225 

10.3 

50 

12967 

1249 

14316 

9.4 

60 

13079 

1216 

14294 

8.5 
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Figure 9 


Effect of Flap Angle on Drag Coefficient 
Change, Referenced to the Drag 
Coefficient for 0° Flap 
Deflection 
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Figure 11 

Experimental XV-15 Download As a Function 
of Flap Deflection Angle 






CHAPTER 3 


Analytic Determination of the Download 

The airflow from a spinning propeller, or prop-rotor 
in the present case, can be approximated through the use of 
vortex ring theory. A vortex ring is a vortex filament 
that forms a circle of radius r'. The approximation of the 
rotor's airflow in this study consisted of a set of concen- 
tric vortex rings, one set for each prop-rotor, beginning 
at a position that coincided with the geometric location of 
the rotor disks. This study was concerned with hovering 
performance simulation only, so the number of vortex rings 
to be included in the analysis was dictated by the hover 
scenario being analyzed, i.e., IGE or OGE. 

The velocity of the air that impinges on the upper 
surface of the aircraft's wing is the cause of the rotor/ 
wing download and hence was the parameter that needed to be 
determined. More specifically, the axial component of the 
rotor slipstream velocity needed to be determined. 

Three aspects of the vortex system that required 
investigation in order for the analysis to be performed 
were: a set of equations that could be used to calculate 

the induced axial velocity components of each of a number 
of vortex rings; an equation that could model the rotor 
wake contraction, i.e., predict the radius (r') of each of 
the vortex rings downstream of the rotor disk; and a system 
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whereby spacing between each successive vortex ring could 
be adjusted to obtain an average induced axial velocity, at 
the rotor disk, that agreed with that predicted by momentum 
theory for helicopter rotor inflow. 


Vortex Ring Induced Velocity 
The velocity induced by a vortex ring can be deter- 
mined by beginning with the Biot-Savart law written in 
vector notation: 


dv 


r ft x ds 
4i ifii 3 


( 10 ) 


giving the velocity vector dv induced by an element of 

length ds = r'd$*' of the vortex ring at the point (x,r,0) 

♦ 

or (x,y,z) (Figure 12). Forming the vector ds/ds , which 
represents the unit vector tangent to the vortex ring at 
the position of the vortex element, with the unit vectors 
in the x, y and z directions gives: 


ds 

ds 


3 


OT ♦ cos#'3 " slnt'ic 


( 11 ) 


Forming the radius vector R, which represents the distance 
from the vortex element to the point (P) at which the 
induced velocity is to be calculated, in a similar fashion 
gives: 


fi * xt + (rslnf - r'sinOj + (rcos# - r'cos**)ic 


(12) 


with the magnitude of the vector being: 
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ifti * j x 2 ♦ r 2 + r' 2 - 


2rr'cos(* -♦') 


(13) 


The cross product of the two vectors then becomes: 


8 * -jf- * -(rcos(* -♦*) - r')T + xs1n#'3 ♦ xcos^’f 


(14 ) 


The x, or axial, component of the induced velocity caused 
by the whole vortex ring can then be determined by inte- 
grating the resulting Biot-Savart expression around the 
circumference of the ring. The expression to be integrated 
becomes : 




V x,r) s 


r 


rcos (»-♦') -1 


[Jx 2 + r 2 + 1 - 2rcosU-*‘) ] 




(15) 


where the ratios x/r 1 and r/r' have been replaced by x and 
r for clarity. Following the steps outlined in Appendix C, 
the final form of the equation to be solved for the axial 
velocity component of the vortex ring is: 


v ix (x.r) 


2ir' 


2 2 
4 x + (r+ir 


K(k) 


1 + 


2(r-l) 

2 2 
x + (r-lT 


E(k) 


( 16 ) 


where K(k) and E(k) are those functions also shown in 
Appendix C. The axial velocity shows singular behavior at 
x = 0, r=r' =1 and near this point the equation can be 
evaluated by expanding the elliptic integrals for K(k) and 
E(k) in powers of (1 - k 2 ): 



m) = A+(A-l)/4(l-k 2 )+(9/64)(A-7/6){l-k 2 ) 2 + 


( 17 ) 


E(k) = l+l/2(A-l/2)(l-k 2 )+(3/16)(A-13/12)(l-k 2 ) 2 + 


(18 ) 


Where : 


A = 



k 2 = 4r 

* 2 2 

x‘ + (r+ir 


The remaining parameter that required attention in 
order to be able to evaluate the induced axial velocity 
equations was the vortex circulation strength, r. 

Beginning with the Kutta-Joukowski law, which relates 
circulation to lift, or in this case thrust, an expression 
for circulation strength can be written: 


Utilizing the dimensionless thrust coefficient from 
standard helicopter terminology, the rotor -produced thrust 
can also be expressed as: 

T ■ •vfo 

( 20 ) 

Equating these two expressions gives: 


IR ■ Vt C t 


( 21 ) 
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or 


rR «= »R 2 V t C t 


( 22 ) 


Differentiating both sides of this expression with respect 
to R gives: 


rdR = 2iRV t C t dR 


(23) 


or 


r • 2.RV t C t 


(24 ) 


This expression, while representing the circulation 
strength about a rotor disk of radius R, has also been 
shown to represent the total circulation strength about the 
axis of rotation due to the trailing vortex system 
(Franklin, 1955). Since this study involved a three-bladed 
prop-rotor system, it was assumed that the circulation 
strength due to each blade could be obtained by dividing 
the total circulation strength, about the axis of rotation, 
by three. Thus, the final form of the trailing vortex 
circulation strength due to each blade becomes: 

2.RV.C. 
r 3 

(25) 


Wake Contraction 

A formula was needed that could be used to model the 


contraction of the rotor wake as a function of downstream 



position from the rotor disk. From the literature investi- 
gated (Marr, Sambell , and Neal, 1973), a formula was 
obtained that, while not reflecting extreme precision, 
modeled the gross characteristics of the wake contraction. 
The expression is of the form 


r(z) 


0.707 


1 


1 1/2 


l 0.5(2 - e' nz/R ) 


(26) 


where: r(z) = downstream wake radius 

R = rotor radius 
z = downstream position 
n = a constant 

The determination of the constant, n, began with the 
continuity equation written between the rotor disk and a 
general point in the rotor wake. It can be shown, between 
these two points, that the following relation holds: 



(27) 


Considering a velocity relationship of the form 



(28) 


it can be written 



(29) 


This relationship was further developed by obtaining 
expressions for the ratios on the right-hand side. For the 
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velocity ratio it is accepted that the ratio of the inflow 
velocity at the rotor disk to that infinitely far down- 
stream is 0.5. Marr, Sambell, and Neal (1973) theorized 
that the radius ratio Rz/Rd equals .908 at a value of z/R 
equal to 0.2. Substituting these values yielded a velocity 
ratio of 

V 

-Y~ = 0.60645 

CD 

(30) 

Considering this relationship and that expressed by Equa- 
tions (26) and (29), it can be shown that 


1/2(2 - e' n2/R ] = 0.60645 
From which it follows that n = 1.19706. 


(31) 


Time Spacing of Vortex Rings 
A method had to be devised whereby the spacing between 
vortex rings could be adjusted so that the average inflow 
velocity through the rotor disk, determined using the com- 
puter simulation model, could be made to agree with that 
predicted by the momentum theory for helicopter rotors 
( i .e. , v. « v -J C fc /2 ) . 

Once again, beginning with the continuity equation, it 
is seen that 


mi . [ R i l 2 

V i " l R(Z) J 


( 32 ) 



Substituting Equation 26 for the radius ratio and squaring 
gives 

J!iiL . ( 2 . ,-■»/*] 

(33) 

Substituting the momentum theory for the induced velocity 
at the rotor disk gives 


V(z) * V t ^7T(2 - e" nz/R ) 

(34 ) 

Realizing that this represents velocity in the z direction, 
i.e., dz/dt, it can be rearranged to the form 


dt = 


1 r dz 

i (2 - *~ n2/R ’ 


(35) 


Integration of this equation yields the final form as 


Kt = 




-f * -§r '<4 ‘ e 


-nz/R 


(36) 


The constant, K, was introduced to provide the means for 
time spacing adjustment of the vortex ring elements, as 
stated earlier. 


Analytic Download Analysis 

With all the components identified and accounted for, 
the task of modeling the XV-15 propulsion system flow, for 
the purpose of analytically approximating the aircraft 
download, could be begun. 



Due to the dependence of the time equation (Equa- 
tion 36) and the velocity equation (Equation 16) itself on 
the rotor tip speed, a constant tip 6peed had to be estab- 
lished. The actual aircraft hovers at 98% of the full 
rotor speed of 601 RPM, or 588.98 RPM. This rotor speed 
corresponds to a tip speed of 768.38 ft/sec and, since it 
represents the actual hover condition, was the constant tip 
speed chosen. At this speed each rotor blade completes one 
revolution approximately every .1022 seconds. The vortex 
system that is actually generated by the rotor resembles 
three continuous shed vortices that together form a con- 
tracting helicoil shape proceeding downstream from the 
rotor disk. It was assumed, for this study, that the 
actual vortex system could be adequately modeled by 
breaking the three concentric helical vortices into three 
concentric vortex rings. In this scheme, three vortex 
rings were considered shed for each rotor revolution, one 
for each rotor blade. It was initially assumed that the 
vortex rings would be equally spaced over the time it took 
the rotor to complete one revolution (i.e., a vortex ring 
was assumed to be shed from the rotor approximately every 
.034 seconds). Because Equation 36 could not be solved 
explicitly for downstream position of the vortex rings as a 
function of time, another procedure was implemented. 

Because the entire analysis was performed for three thrust 
coefficients. Equation 36 was plotted using each thrust 
coefficient and whole values of z, the downstream position 
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of the vortex rings (Figure 13). Then, using five sets of 
points from each resulting curve, a fifth order inter- 
polating polynomial was formed for each thrust coefficient 
dependent curve. This resulted in three fifth order 
polynomials that were explicitly solvable for the vortex 
ring downstream position as a function of ring time 
spacing. Up to this point the time spacing factor, K, had 
been held constant at unity and was retained in the inter- 
polating polynomials. Then, beginning with the initial 
time spacing of the vortex rings (.034 seconds between 
rings), the downstream distance that each ring had been 
displaced from the rotor disk could be calculated. With 
this distance known, the associated vortex ring radii could 
be calculated (Equation 26) and the contribution of each 
ring to the total axial velocity at the points of consid- 
eration could be computed. The ring spacing factor K could 
then be adjusted and the computations performed again until 
successive iterations arrived at the desired result. 




calculations were begun 2.5 feet outboard of the rotor tip 
path and proceeded in one-foot increments through the rotor 
plane to a point 2.5 feet outboard of the rotor tip path on 
the opposite side. For each point the axial velocity 
component due to each of the vortex rings was calculated. 
These components were then summed to arrive at a total 
axial velocity for each point. Because the aircraft rotors 
are 25 feet in diameter, this scheme led to a total axial 
velocity value for each of 31 points. The 25 velocity 
values that occurred in the "rotor disk" were then averaged 
to arrive at an average inflow velocity through the rotor 
disk. The vortex ring spacing factor K was then varied 
until the average inflow velocity through the rotor disk 
agreed with that predicted by momentum theory. Small K 
values packed the rings more closely and resulted in higher 
average velocities, while large K values spaced the rings 
out and resulted in lower average velocities. The range of 
K was between zero and one. The K values that resulted in 
the correct average inflow velocities were then retained 
and used later in the version of the program that computed 
aircraft download. 

The number of vortex rings in the calibration program 
analysis was limited to 72. This number of rings produced 
adequate results evidenced by the fact that increasing the 
number of rings included in the analysis to 100 increased 
the average inflow velocity through the rotor disk by a 
maximum of .39%. 



Download Computer Program 

With the proper predetermined spacing factor K, the 
program version that calculated the analytic download could 
be formulated. 

The geometry of the analysis consisted of two systems 
of vortex rings, each concentric, their centerlines sepa- 
rated by approximately 32 feet. The axial velocities were 
calculated at 71 points; beginning at a point approximately 
6.5 feet outboard of one rotor tip path and progressing, in 
one-foot increments, to a point approximately 6.5 feet 
outboard of the other rotor tip path. The points at which 
the velocities were calculated were 4.67 feet below the 
rotor planes, coinciding with the position of the wing on 
the aircraft itself. 

Both IGE and OGE hover scenarios were investigated. 

For the OGE scenarios (H/D = infinity) , each rotor was 
modeled by 72 concentric vortex rings for a total of 144 
rings. For the IGE scenarios (H/D = 2.0, 1.0 and 0.5), the 
number of vortex rings used was dictated by the H/D ratio. 
For these cases the corresponding rotor heights were 50, 25 
and 12.5 feet, respectively. The ground plane was modeled 
by using a mirror image of the shed vortex system, 
reflected about the "ground plane." This was accomplished 
in the following manner. The H/D ratio to be investigated 
was a program -prompted input. If it was an IGE hover 
condition, this input set the required distance that sepa- 
rated the rotor plane from the ground plane. Then, while 



going through the program loop that calculated the vortex 
rings' downstream positions, the positions were monitored 
until a position was reached that corresponded to the 
required distance between the rotor and ground plane. At 
this point, the next vortex ring that was formed was iden- 
tical to the one that occurred at the ground plane position 
and in fact was superimposed on the ring that occurred at 
the ground plane. The next ring that was formed was iden- 
tical to the ring that preceded the ring that occurred 
before the ground plane position. Each successive ring 
that was formed for the reflected ring system was identical 
to its corresponding ring in the shed vortex system 
(Figure 14). The circulation strengths on the superimposed 
mirror image system, however, were reversed. With all the 
positions and radii of the vortex rings determined, the 
total axial velocity at each point of concern could be 
computed. The velocities from successive points were then 
averaged and used to form an average dynamic pressure, q, 
based on standard sea level conditions. 

Values for the wing drag coefficients were taken from 
wind tunnel tests of a representative, two-dimensional wing 
section. The wing section was equipped with movable flaps 
that permitted testing at deflections of 0°, 20°, 40°, 45°, 
and 60°. The wing section was tested at angles of attack 
near -90°. The results of these tests provided five differ 
ent drag coefficients corresponding to the five flap angle 



43 



deflections. The flap deflection angle to be analyzed was 
a prog ram -prompted input. 

With the averaged dynamic pressures and drag coeffi- 
cients determined, all that remained to be done, for 
download or drag considerations, was to determine an appro- 
priate surface area. Due to the fact that the points used 
to arrive at the averaged dynamic pressure were separated 
by one foot, this one-foot length was multiplied by 5.25 
feet, the aircraft wing chord. This 5 . 25-6quare-foot area 
was used between the vortex system centerlines, as the 
rotors are positioned on the wing tips of the aircraft. 

Each averaged dynamic pressure was then multiplied by a 
5. 25 -square -foot -area and the appropriate drag coefficient, 
depending on the flap deflection angle, to arrive at a 



component of drag or download. The drag components across 
the wing were then summed to arrive at a total value of 
drag. The resulting total drag or download value was then 
divided by the total rotor -produced thrust, determined by 
the thrust coefficient being investigated. This gave 
results for download as a percentage of total aircraft 
thrust. Listings of both the calibration and download 
versions of the XV-15 program are included in Appendixes D 
and E. 


Velocity Profiles 

During the formulation of the propulsion system com- 
puter simulation program, it was necessary to check each of 
the velocity components against published values (Deckert 
and Ferry, 1960) before the total velocities were computed. 
This was done in order to check the accuracy of the values 
from the computer program. The published values are 
presented as numbers that are dependent on the rotor radius 
to vortex ring radius (r/r 1 ) ratio and also on the distance 
between the point at which the velocity is being calculated 
and the vortex ring centerline to vortex ring radius (x/r') 
ratio. The published numbers are then multiplied by two 
pi, the vortex ring radius ( r ' ) and divided by the circula- 
tion strength to arrive at the velocity value. Manually 
calculating these two ratios and knowing the circulation 
strength permitted a spot check of some of the number of 
velocity components that resulted from the computer program 
against the published values. When these values agreed, it 



was concluded that the program was functioning correctly 
and the correct total velocity at each point would result. 
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An interesting intermediate result was obtained when 
the program was next set up to calculate the total axial 
velocity at each point along the line between the two rotor 
centerlines. These results produced a two-dimensional 
cross section of what was theoretically occurring to the 
rotor airflow due to the proximity of the rotors and the 
geometry of the installation. The resultant total axial 
velocities were bilaterally symmetric and are presented for 
each of the three thrust coefficients in Figures 15, 16, 


and 17 . 
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Figure 17 

Computer -Generated Velocity Profile, 
C t = .0105 


CHAPTER 4 


Model Test 

After it was determined that the velocity values from 
the computer program were accurate, it was desirable to 
formulate a test for the download aspect of the program to 
demonstrate its validity. This was accomplished by using a 
model rotor system, a schematic of which is presented in 
Figure 18 . 



Figure 18 

Model Rotor Test Setup 

The model consisted of a 1.135-foot radius, fixed 
collective pitch rotor. The rotor had no provisions for 
cyclic pitch. The rotor was powered by a small DC motor 
that was powered by an instrumented power source. Power 
was transferred from the motor to the rotor by a shaft and 
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gear system. The power that was input to the motor was 
monitored through the use of volt and ampere meters on the 
power source. The thrust output by the rotor was deter- 
mined through the use of a spring scale mounted at the 
opposite end of a statically balanced lever arm. 

The model rotor system was operable both with and 
without a flat plate "wing" section placed 1.5 inches below 
the rotor disk. It was thought that by operating the model 
system both with and without the wing installed, monitoring 
the power source, spring balance and rotor RPM, the down- 
load on the flat plate "wing" could be deduced. 

A computer program was formulated to simulate the 

f 

model performance. Both the calibration and download 
versions of the model program are included in Appendixes F 
and G. These programs were constructed in exactly the same 
manner as were the programs that modeled the XV-15. 



CHAPTER 5 


Results and Discussion 

Experimental Download Analysis 
Values of rotor/wing download for the XV-15 Tilt Rotor 
Research Aircraft were deduced by combining the results 
from actual aircraft hover tests and wind tunnel tests of a 
representative two-dimensional wing section. 

The hover tests were conducted at NASA Ames Research 
Center in 1983. The aircraft was modified to allow the 
flap and flaperon control surfaces to operate as one. This 
provided the effect of a full span plain flap system. The 
aircraft gross weight was closely monitored and maintained 
within a 200 -pound range. The rotor speed was held con- 
stant during the hover test points. 

The wind tunnel tests were also conducted at Ames 
Research Center, but in the Army 7-foot by 10-foot wind 
tunnel. These tests resulted in drag coefficients for an 
airfoil section that had a .25c operable flap and was 
tested at angles of attack near -90°. 

The values of the rotor/wing download that resulted 
from combining the data of these two tests are tabulated in 
Table 5 (page 24) and presented graphically in Figure 11 
(page 27). The values determined for the download phenom- 
enon ranged from 13.2% of total aircraft thrust (DL/T) for 
a 0° flap deflection to 8.5% (DL/T) for a 60° flap 
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deflection, both at a power coefficient (Cp) equal to .001. 
While these results appear reasonable, they are larger than 
those previously predicted for the XV-15 download phenom- 
enon. Investigation of past studies that analyzed the 
download phenomenon, however, revealed that the magnitude 
of the resulting download values has increased. Thus the 
larger values produced by this study agree with the trend. 
However, there are some aspects of the experimental anal- 
ysis that merit investigation. 

The XV-15 was not flown in its regular configuration 
for the hover tests used for this study. A special modifi- 
cation was made to the flap/flaperon system so they could 
be deflected as one control surface. Thus, the data taken 
for this study reflect an XV-15 aircraft with a full span 
plain flap instead of the actual flap/flaperon system. 

While some data exist that reflect the actual flap/flaperon 
system, the limited number of test points and data scatter 
prevent a comparison to the full span plain flap data. 

Hence, as a result of this study, it can be deduced that 
the actual download for the XV-15, in the nominal hover 
configuration (40° flap/25°f laperon ) lies between the 25° full 
span plain flap value of 11.5% DL/T and the 40° full span 
plain flap value of 10.4% DL/T. Hence, more OGE hover 
tests of the XV-15 with the "as-built" flap/flaperon sys- 
tem, at various angles of deflection, are required to 
narrow this download range and better define the XV-15 


download . 



The wind tunnel tests of the two-dimensional wing 
section, while providing valid results for the "as tested" 
configuration, did not adequately model the actual flow 
field in the vicinity of the aircraft. Since the tests 
have been completed, more understanding has been developed 
about the flow field that exits about the XV-15. A region 
exists, midway between the rotor centerlines, where there 
actually is an "upwash." This is caused by the radial 
components of the rotor wakes that travel along the wing 
span and intersect in the vicinity of the aircraft center- 
line, turning upward. This flow phenomenon was validated 
in the subsequent computer simulation model of the XV-15 
done for this study. The existence of this upwash is 
documented in the velocity profile plots (Figures 15, 16, 
and 17). Due to this fact, and another more obvious limi- 
tation, the extension of the two-dimensional drag data to 
adequately represent the aircraft wing drag is question- 
able. The more obvious limitation is encountered in the 
process of transforming two-dimensional airfoil section 
data to reflect three-dimensional characteristics. While 
the transformation of sectional data to reflect wings of 
finite aspect ratio is, in general, easily accomplished, it 
is not so in this case. A better representation of the 
airflow about the wing of the XV-15 must be established in 
order to obtain more closely representative drag data. It 
is believed that tests must be performed with a wing, 
instrumented to obtain force and moment data, placed 



downstream and in between two contra-rotating prop-rotors. 
Whether this test is run with scale-model or actual-size 
hardware, it would further the understanding of the wing 
drag characteristics as they apply to the XV-15 configura- 
tion. 

The assumption that the download-to-thrust (DL/T) 
ratio was a constant for a given flap angle deflection was 
found to be accurate when applied to a constant Cp anal- 
ysis. The download-to-thrust ratio was found to actually 
increase slightly with increasing power coefficient, but 
this increase amounted to only 2% over a range of Cp 
from .0006 to .0015. 

A method that more accurately predicts the engine 
exhaust jet thrust must be developed. The method used for 
this study relied on the assumption that the average inflow 
velocity, as predicted by simple momentum theory, ade- 
quately represented the velocity in the vicinity of the 
engine inlet. To better understand the engine inflow 
velocity, pressure surveys of that area should be per- 
formed. The magnitude of the exhaust thrust, however, was 
small enough that its contribution to the analysis was 
mi nor . 

3/2 

The correction of the power polar (Cp vs . C fc ) 
algorithm to reflect the nonlinear aerodynamic effects at 
high Cp values also deserves more study and validation. 



Model Test Analysis 


The model test and subsequent computer simulation of 
the model did not supply results with the desired preci- 
sion. The results that were produced, however, compare 
favorably. The model test resulted in a download-to-thrust 
(DL/T) ratio of 11.2% of rotor thrust at a thrust coeffi- 
cient of .0075 and an H/D ratio of infinity. The computer 
simulation of the model, however, produced a DL/T ratio of 
approximately 9.7% for the same conditions. The difference 
(approximately 13%) between the two values could be attrib- 
utable to a few factors. 

First is the validity of taking two-dimensional drag 
data and applying it to a three-dimensional test. The drag 
coefficient for the computer simulation was taken to be 
1.98 based on the literature reviewed (McVeigh, 1985). 

This represented a flat plate oriented perpendicular to the 
flow (i.e., an angle of attack of -90°), tested at Reynold' 
numbers between 10,000 and 1,000,000. Since there were no 
provisions on the model setup to measure the velocities 
that occurred at the wing surface, a determination of the 
test Reynold's number was not possible. It is also thought 
that a better understanding of the drag characteristics 
should be obtained. This would allow the computer program 
to model more precisely the actual test conditions. A 
model system should be built and tested that incorporates a 
wing, instrumented to provide force and moment data, sepa- 
rate from the rotor force system. This would allow direct 



measurements and correlation of the download force on the 
wing and the rotor-produced thrust. 

While some of the aspects of the model test could 
benefit from further refinement, the test provided an 
interesting result. The comparison between the model test 
and computer simulation did show, in this case, that the 
computer simulation provided a lower estimate of the down- 
load than did the model test. This provided the valuable 
information that the computer simulation predicts lower 
download values than those which actually occur. Extended 
to the aircraft simulation program, this means that the 
magnitude of the actual download is at least as large as 
the analytically determined download. 

Analytic Download Analysis 

This aspect of the study was developed with contribu- 
tions from a number of sources. The main component of the 
study was the theory involved in modeling an aircraft rotor 
system's airflow through the use of vortex ring theory. 

Formulas were developed that modeled the rotor wake 
contraction downstream of the rotor disk (Equation 26) and 
also allowed for adjustment of the time spacing of the 
vortex rings in the wake (Equation 36). A formula was also 
obtained that allowed determination of the needed vortex 
ring circulation strength (Equation 25). Then, beginning 
with the Biot-Savart law, a formula was developed that 
permitted calculation of the axial velocity component due 
to each of the vortex rings (Equation 36). The velocity 



components were then summed and used in conjunction with 
the aforementioned wind tunnel drag coefficients to arrive 
at download values. This analysis was performed for three 
different thrust coefficients and four H/D ratios for each 
thrust coefficient. These download values are tabulated in 
Table 6. Figure 19 is presented to show the effect that 
different H/D ratios had at a constant thrust coefficient 
of .0094 on the theoretical download values. The effects 
that different thrust coefficients had on the download 
results were too small to allow an adequate graphic com- 
parison. The range of download values for this method of 
analysis varied from 14.44% DL/T for a 0° flap deflection to 
9.233% DL/T for a 60° flap deflection, both at H/D equal to 
infinity and equal to .0105. Changing the H/D ratios 
from infinity to 0.5 had the effect of reducing the DL/T 
ratio by from 5% to 7%. There are, however, a few aspects 
of this method of analysis that warrant further investiga- 
tion . 

The method by which the circulation strength of the 
vortex rings was determined requires further investigation 
and refinement. The velocity components that were computed 
using this method are directly proportional to the circula- 
tion strength, hence the circulation strength is an 
important parameter. While the method presented seems 
rational, it is thought that other factors that were not 
investigated could affect the circulation strength. One 
such aspect is the assumption that the circulation strength 
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Table 6 

Analytic Download 


c t 

h/D 

tI **> 
6 f =0° 

tI <*> 

1 J f = 20 ° 

tI <*' 

« f s 40 ° 

tI <*> 

« f = 45 ° 

rl <*> 
"« f = 60 c 


00 

14.429 

12.919 

11.241 

10.737 

9.227 

.0094 

2.0 

13.848 

12.403 

10.797 

10.315 

8.869 


1.0 

11.919 

10.691 

9.325 

8.916 

7.687 


0.5 

7.368 

6.632 

5.814 

5.569 

4.833 


00 

14.424 

12.914 

11.236 

10.733 

9.223 

• 

o 

H* 

2.0 

13.807 

12.366 

10.765 

10 . 285 

8.844 


1.0 

11.717 

10.511 

9.170 

8.768 

7.562 


0.5 

7.567 

6.811 

5.970 

5.718 

4.961 


00 

14.440 

12.928 

11.248 

10.744 

9.233 

. 0105 

2.0 

13.784 

12.345 

10.747 

10.268 

8.829 


1.0 

11.794 

10.579 

9.229 

8.824 

7.609 


0.5 

6.625 

5.967 

5.236 

5.017 

4.360 






remains constant in the vortex ring as it translates down 
the rotor wake, while it is believed that a decay of circu- 
lation strength actually occurs. In order to better 
resolve this issue, rotor tests that address the circula- 
tion strength phenomenon should be performed. 

The formula (Equation 26) that models the rotor's wake 
contraction as a function of downstream distance from the 
rotor disk also requires more investigation. Analytical 
methods have been produced since Equation 26 was developed 
in 1982 that more closely model the actual behavior of the 
wake downstream of the rotor disk. These methods should be 
incorporated into this type of computer simulation analysis 
to reflect the increased knowledge of rotor wake 
contraction. The constant, n, in Equation 26 should be 
investigated further to determine its control of the final 
wake shape and overall characteristics. It is thought that 

* 

a better method of solving for this constant could be 
devised and implemented in the analysis. 

The fact that this computer simulation concerned 
itself with only the airflow due to the two aircraft rotors 
and not any of the contributions due to the body effects of 
the aircraft itself sheds some doubt on the precision of 
the final answers. While the resulting values of download 
are reasonable, the inclusion of the aircraft body into the 
analysis would no doubt have an effect on the final 
results. For further study, the results of this analysis 
could be used in conjunction with a method that included 



the physical body of the aircraft to better understand the 
wake/body characteristics that contribute to the download 
phenomenon. This could also be extended to gain a better 
understanding of the sensitivity of the download to basic 
configuration changes of the aircraft body shape. 

While extension of the wind tunnel data to reflect the 
actual aircraft conditions is suspect, in the case of the 
computer simulation it is thought that the wind tunnel 
results adequately represent the situation that exists in 
the computer simulation. The computer simulation concerned 
itself with only the determination of the axial components 
of the velocity that exist on the upper wing surface. 

Since the wind tunnel tests essentially reflected only 
axial and not radial velocities, i.e., spanwise flow, it is 
thought that the results from the wind tunnel tests and 
computer simulation were compatible. Hence, the results 
from both tests, when combined, were thought to form a 
valid analysis procedure. Further tests should be per- 
formed to either validate this assumption or develop a more 
nearly correct means of representation of the flow charac- 
teristics . 

The intermediate results of the velocity profiles 
demonstrated the validity of the upwash phenomenon that 
occurs between the rotors along the aircraft centerline. 
This upwash has since been further validated through the 
use of flow visualization techniques applied to both the 
aircraft itself and model tests. The inclusion of this 



phenomenon in the download picture presents another factor 
to be considered. 



CHAPTER 6 


Conclusion 

Analysis of the download phenomenon, which occurs when 
the XV-15 Tilt Rotor Research Aircraft is flying in the 
hover/helicopter configuration, was pursued utilizing two 
different methods. The download is a result of the air- 
craft's wing being in close proximity (approximately 5 
feet) and parallel to the two rotor disks. The download 
force has the direct effect of reducing the aircraft's 
useful hover load and hence is an important parameter when 
considering tilt rotor aircraft. 

The first analysis method utilized results from actual 
aircraft hover tests and wind tunnel tests of an XV-15 
representative wing section. The download results for this 
method ranged from 13.2% DL/T for a 0° flap deflection to 
8.5% DL/T for a 60° flap deflection, both at an H/D of 

infinity and a C of .001. 

P 

The second analysis method relied on an analytical 
model of the aircraft's two prop-rotors and their asso- 
ciated airflow. The model was then used to determine the 
axial velocities that occur at the position which the 
aircraft wing occupies. By then applying the drag coeffi- 
cients from the aforementioned wind tunnel tests, download 
values resulted. A computer program was utilized for this 
analysis method. The results for this analysis method 
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ranged from 14.44% DL/T for a 0° flap deflection to 9.233% 
DL/T for a 60° flap deflection, both at an H/D of infinity 
and a C t equal to .0105. The analytical analysis was 
performed for each of three thrust coefficients (.0094, 

.01, .0105) and also included the investigation of four 
H/D (infinity, 2.0, 1.0, 0.5) ratios for each thrust 
coefficient. The effect that the different H/D ratios had 
on the download results was quite apparent. Reducing the 
H/D ratio from infinity to 0.5 had the effect of reducing 
the DL/T ratio by from 5% to 7%. However, the effect that 
the different thrust coefficients had was not pronounced 
enough to allow any conclusion to be drawn. 

A model rotor system was tested and analytically 
modeled in an attempt to validate the computer analysis. 
While the results did not agree with the desired precision, 
reflecting a 13% difference, the download results of the 
computer analysis were lower than those of the model test. 
This implied that the computer program predictions were 
within the lower bounds of the actual aircraft download 
values . 

This study does not represent the definitive download 
analysis. At the time that this report was written, there 
were at least two other separate studies being conducted in 
attempts to better define the XV-15 download issue. One 
study is utilizing a powered model of the XV-15 aircraft 
with wing surfaces instrumented to obtain force and moment 
data. An additional objective of this test is to 



investigate flow control, i.e., suction and blowing, and 
its ability to lessen the model download forces. Another 
study, utilizing the XV-15 aircraft, is attempting to 
establish a relationship between rotor thrust and hub 
spindle beam bending moments to resolve the download issue. 
It is thought that if this relationship can be established 
from data taken with the aircraft on the tie-down stand at 
NASA Ames Research Center, it can be applied to data 
recorded in free hover flight and download values will be 
determinable. 

The download phenomenon is an extremely important 
performance parameter for the tilt rotor family of air- 
craft. When viewed in a useful load context, the download 
can be responsible for large reductions in the aircraft 
hover weight capabilities. With the V-22 aircraft coming 
into production for military and even possibly civil appli- 
cations, more effort must be made to better understand and 
define the download phenomenon. It is only when a problem 
is fully understood and defined that intelligent attempts 
to deal with it are possible. 
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Transformation from Equation 6 to Equation 7 
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Appendix B 


Empirical Adjustment to Power Polar to 
Reflect Nonlinear Aerodynamics Region 
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Appendix C 

Development of the Biot-Savart Equation 
to Determine the Axial Velocity 
Component Due to a Vortex Ring 
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Appendix D 

Calibration Version of XV-15 
Computer Program 
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ORIGINAL PAGE IS 
OF POOR QUALITY 




C 

C 

C THIS PROGRAM IS DESIGNED TO SIMULATE THE XV- 15 ROTOR SYSTDTlTERFORMARCE 
C THROUGH THE USE OF A VORTEX RINC ANALYSIS. THIS PROGRAM IS! BASED ON 

C A CONSTANT ROTOR TIP SPEED OFl 768.38 FT/SEC AND CAN PERFORM VTHE 

C ANALYSIS FOR THREE DIFFERENT THRUST COEFFICIENTS ( . 6094 ,.611.91 03 ) . 

C THIS PARTICULAR PROGRAM VERSION WAS DESIGNED TO CALIBRATE THE VORTEX 

C RINC SPACINC FACTOR • (K> AND HENCE MODELS ONLY ONE ROTOR OF THE SIDE 

C BY SIDE ROTOR SYSTEM. THE VALUE FOR IC IS AN INPUT PARAMETER IN THIS 

C PROGRAM AND IS VARIED UNTIL AN AVERAGE INFLOW VELOCITY. THAT AGREES 

C V/77? THAT PREDICTED BY 8IHPLE MOMENTUM THEORY. V*VMp!Ct/2. >**.3. 

C IS OBTAINED. 

C DUE TO THE FACT THAT THE AVERACE INFLOW VELOCITY AT THE ROTOR DISK 
C IS WRAT*S DESIRED. THE POINTS 1 AT WHICH THE VELOCITIES ARE CALCULATED 

C LIE ALONG A LINE WHICH IS COINCIDENTAL WITH THE ROTOR PLANE. 

C 

C 

M *ammm*MM* MMm*mm MmMM*9 mammm*9*m*******m**mM*99*Mma Mm**m*» ** 


r**»»»***» tnsmttsmu 


«*anr**tx**xx*****«xKaa 




LIST OF VARIABLES: 


K VORTEX RINC SPACINC FACTOR 

B ROTOR SPEED ! SEC/ REV) 

F ( I > TIME BETWEEN VORTEX RINGS fSEC) 

(BASED ON THREE RINGS PER REVOLUTION) 

KT(I) CORRECTED TIME BETWEEN VORTEX RINGS (SEC) I 

(K*F( I ) ) 

Z< I > THE POSITION OF EACH VORTEX RING. DOWNSTREAM 

OF THE ROTOR DISK. AS A FUNCTION OF CORRECTED 
TIME (FT) 

RR(I) VORTEX RING RADIUS (FT) 

(A FUNCTION OF ROTOR GEOMETRY AND DOWNSTREAM 
POSITION) 

TNR. . . THE TOTAL NUMBER OT VORTEX RINGS CONSIDERED 

IN THE ANALYSIS 

A THE VORTEX RINC CIRCULATION (FT**2/SEC> 

CT THRUST COEFFICIENT ( THRUST/ ( RHO* AREA* V t t p**2 . >> 

L THE DISTANCE BETWEEN THE ROTOR AXIS AND THE 

POINT WHERE THE AXIAL VELOCITIES ARE BE1NC 
CALCULATED (FT) 

J THE ABSOLUTE VALUE OF POINT (FT) 

X THE RATIO OF VORTEX RINC DOWNSTREAM POSITION 

TO RING RADIUS (Z(I>/RR(D) 

R THE RATIO OF J TO THE: RING RADIUS 

(RAD/RR( I ) ) 


9. T.U.E.B.C. CIRC. .COMPONENTS OF THE EQUATION USED TO DETERMINE 


THE AXIAL VELOCITY DUE TO EACH RINC ( FROM: 
'AERODYNAMICS OF PROPULSION’ .KUCHEHANN 8 WEBER. 
1953) 

P.D INTERMEDIATE STORACE VARIABLES 

V( M , H ) THE COMPONENT OF THE AXIAL VELOCITY AT POINT 

H DUE TO RINC N (FT/SEC) 

VTOT(l) THE TOTAL AXIAL VELOCITY AT EACH OF 31 POINTS 

DUE TO SUKMINC THE COMPONENTS FROM EACH VORTEX 
RING (FT/SEC) 
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C 

C 

C 

c 

c 

c 


VAVC. . 

vdist. 


c 

c 


THE SOT! Of THE AXIAL VELOCITIES'. ALONG THE ROTOR 
disk < nrvsEc > 

THE AVERAGE: INFLOW VELOCITY ALONG THE ROTOR DISK 
< FT/SEC ) 


** ** * ** *** **** 


INTECER TNR 

REAL K.KT<72> ,RR!72> ,Z!72) ,FC72) ,VT0T!3l ) ,V<3l ,72) 
C 
C 






************* 


C 

C 

C THE VALUE TO BE USED FOR THE VORTEX RING SPACINC FACTOR CK) IS 
C FIRST READ IN AND THE USER IS PRQHPTED FOR THE DESIRED THRUST 
C COEFFICIENT <CT>. 

C 

C 

C ******** a* »»a**y«a* * g«»*a* «*«a** **»a***«*«***s« « »*« ««*** 

C 

C 

VRITE<6,4> 

4 FORMAT! * WHAT K VALUE WILL BE USED?* ) 

R£AD< S.3)K 
3 FORNAT ( C 1 2 . 5 > 

WRITE(B.9) ' 

9 FORMAT!’ WHAT THRUST COEFFICIENT (CT) WILL BE USED?* ) 

WRITE! 6 *11) 

11 FORMAT! * ( .0094, .01 , .0103)?* ) 

READ! 5, 10 )CT 

10 FORMAT ! Cl 2 . 3 ) 

C 

C 

C 

C 

C NEXT. BASED ON A ROTOR SPEED OF 760.38 FT/SEC !.1022 SEC/REV) . THE 

C TIME SPACING OF TUEl VORTEX RINGS IS CALCULATED USING THE INPUT RINC 

C SPACINC FACTOR !K). 

C 

C 

Cttmta»x»>uu»u»»»>iu<xxitinnx(mtxtt»xixx<tnMSX«*uxtxxiLxxxtt 

C 

C 

H= . 1022 
DO 20 1=1.72 

F ( I > « I aH/3 . 0 
KT ! I ) ■ K* F < I ) 

20 CONTINUE 


C 

C*i 

C 

C 

C 

C 

C 

C 


ixxaxxxxmxxi 


***************** 


DEPEND INC ON WHICH THRUST COEFFICIENT WAS SELECTED. AN IF/TBEN 
STATEMENT SETS THE PROPER VALUE OF VORTEX RINC CIRCULATION AND 
TRANSFERS CONTROL TO THE APPROPRIATE INTERPOLATING POLYNOMIAL. 
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C 

G 

(F(CT.Ea..e*94> THEN 
A* I 89. 08 
CO TO 25 
ENDIF 

IF(CT.EQ. .61 ) THEN 
A*201 .10 
CO TO 39 
END1F 

IFtCT.EQ. .0163) THEN 
A*2l I . 22 
CO TO 33 
END IF 


C 

C 

c 

c 

C D EPEN D INC ON WHICH THRUST COEFFICIENT WAS SELECTED, ONE OF THREE 

C INTERPOLATING POLYNOMIALS 19 USED TO CALCULATE THE DOWNSTREAM 

C POSITION OF EACH VORTEX RING AS A FUNCTION OF TIKE. EACH VORTEX 

C RING'S RADIUS IS ALSO CALCULATED AS A FUNCTION OF RING DOWNSTREAM 
C P09ITI0N . 

C 

C 

c*»********** ******************»********«*************** **** ***** * ********** 

c 

c 

23 DO 46 1*1.72 

Z< I ) »-63 . 833*KT< I )**3 . *200 . 424*KT< I >**4 . -242 . 399*ICT< I >**3 . 

♦ ♦ 1 20 . 798*KT ( I >**2 . ♦75. 9 12*KT( I > 

RR( I ) * 8 . 837 3*SQRT ( I .✓< . 5* <2 . -EXP<- .09312*24 1 ) ) > ) ) 

40 CONTINUE 

CO TO 24 
C 
C 
C 

39 DO 43 M .72 

Z < I )*-6B. 152*KT< I >**3.4221 .77*KT< I >**4 . -236 .836*KT( I >**3. 

♦ 4126.62*KT( l >**2.470.083*ICTU ) 

RR< I ) * 8 . 837 3*SQRT ( I . 3*< 2 . -EXP ( - .093 l 2*Z< !>>>)> 

49 CONTINUE 

CO TO 24 
C 
C 
C 

39 DO 60 1*1.72 

Z< I )»-64.746*KT< I >**5.42! I .736*KT( I >**4. -247 . 2 1 9*ICT( I >**3. 

4 4123. 687*KT ( I >**2 . +81 .398*KT( I ) 

RR( I > *8. 837 3* SORT ( 1 .✓< . 3* (2 . -EXP ( - .09312*Z( I > > > ) ) 

30 CONTINUE 

C 

C 

C************* **************************************************** ************ 

C 

c 

C THE TOTAL NUMBER OF RINGS TO BE USED FOR THE ANALYSIS (TNR) AND 

C THE POSITION OF THE I FIRST POINT. AT WHICH THE AXIAL VELOCITIES ARE 
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C TO BE CALCULATED. 

C 

C 


c 

c 

24 TMR-72 
L*-l 3 


WITH RESPECT TO TOE ROTOR AXIS ARE READ IN. 


C 

C 

C 

C 

C TWO DO LOOPS ARE USED TO CALCULATE EACH VORTEX RING’S AXIAL 
C VELOCITY COMPONENT AT EACH POINT OF INTEREST. THE INNER LOOP 

C CALCULATES THE VELOCITIES FROM EACH OF TOE RINCS AND TOE OUTER LOOP 

C INCREMENTS THE POINT AT VHICHiTBE VELOCITIES ARE CALCULATED*. 

C 

C 

m tmm mmt iKtwtfWMMmMmtwMWw 

C 

C 

DO 63 H* I .31 

C 

DO 70 N* 1 .TNR 

IF(L.LT.0)J*~L 

1F<L.CE.6U*L 

X«-Z(N>^RR<N> 

R*J/RR(N) 

8«<4.*R)/<X**2.<HA*I . )**2. ) 

T*AL0C<4.^(1 .-S>**.3> 

U«Tf((T-l . )/4. )«< 1 . —S ) ♦ ( 4 . /'64 . )* IT- <7. ✓«.)>*( I -S>**2. 

E*1 .♦< . 3* ( T— . 3 1 * ( 1 ,-S) )M3./16. )*<T-< 13./I2. > >*<l .-S>**2. 

D* I ./<X**2.+<R+1 . >**2. ) ** . 3 

C*-< I .♦( I2.*(R-I . ) >/'<X**2 . ♦<R-1 . >**2. ) ) ) 

D*U*C*E 

P*B*D 

Cl RC 3 A/' < 2 . *3 . 1 4 1 59*RR < N ) ) 

V(H.N)-CIRC*P 
70 CONTINUE 

C 

L=L* 1 

C 

63 CONTINUE 

C 

C 

Ot»«tu*u»»atsnnu*»sia«uttxsnm«<«sx>nnt*suxxiaRm*a*<*«aoitxx» 


C 

C 

c 

c 

c 

c 

c 

C 1 


THE VELOCITY COMPONENTS DUE TO EACH OF THE VORTEX RINCS ARE NO L * 
SUMMED TO ARRIVE AT A TOTAL AXIAL VELOCITY AT EACH POINT OF 
INTEREST 


tmtanytumtutMM 


isms 


C 

C 

DO 73 J* I .31 

DO 76 1*1 .TNR 

VT0T< J )*VT0T( J )+V< J. I ) 
76 CONTINUE 
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73 CONTINUE 

C 

C 


C 

c 

C THE TOTAL AXIAL VELOCITY AT EACH POINT OF INTEREST IS ROW PRINTED 

C OUT ALONG WITH THE ASSOCIATED VALUE OF VORTEX RING SPACING FACTOR 

C (K) . THRUST COEFFICIENT CCT) , AMD THE HUMBER OF VORTEX R1NCS USED. 

C 
C 

C 

C 

VRITE< 6 .83 )K.CT.TNR 

85 FORMAT! IX. * FOR K** ,F7. 6. * « CT*’.F3.4. 

♦ * AND* . IX. 12 . * VORTEX RJNCS* ) 

VRITEI6 90) 

96 FORMAT! IX. * THE TOTAL AXIAL VELOCITY AT EACH POINT ALONG A LINE* ) 
WRITE! 6 . 95) 

95 FORMAT! IX.* COINCIDENTAL WITH THE PLANE OF TOE ROTOR IS * * ) 


DO 97 J* I .31 

WRITE! 6 . 96 ) J . VTOT! J ) 

96 FORMAT !/,2X. 3 HVTOT! ,12.2U)= .C12.3) 

97 CONTINUE 
C 

C 

. C**********m*M*n*********.m**9*mm**9m***nm*Mm**9ic*******m*m***9 ***Mx**wm*rwar* 9** 

C 

C 

C SINCE FOR CALI ORATION PURPOSES THE AVERACE VELOCITY ALONC THE I\OTOn 
C DISK IS WHAT’S OF PRIME IMPORTANCE. THAT VALUE IS NOW CALCULATED 

C AND OUTPUT. 

C 

C 

C«xt«um»m*umxt«aiwmtm«mixx«nstssn»xmnxtntt*ix«xs«mx>u 

C 

C 

DO 69 1*4 ,28 

VAVC* VAVC+ VTOT ( I ) 

69 CONTINUE 

VDlSK=VAVC/23.0 
VRITE(6.7I )VDISK 

71 FORMAT!/.’ THE AVERACE I NTLOW VELOCITY IS* . IX.F5.2. * FT/SEC') 

STOP 

END 



Appendix E 

Download Version of XV-15 
Computer Program 
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C 

C 

C THIS PROCRAN 19 D ESI CITED- TO SIMULATE THE PROPULSION SYSTEM OF THE 

C XV- 13 AIRCRAFT. THROUGH THE USE OF A VORTEX RING ANALYSIS. FOR THE 

C PURPOSE OF DETERMINING THE ROTOR GENERATED WING DOWNLOAD THAT OCCURRS 

C IN THE HOVER CONFIGURATION. THE PROCRAM IS BASED ON A CONSTANT ROTOR 

C TIP SPEED OF 768.30 FTVSEC AND UTILIZES' THE PREDETERMINED VORTEX RING 

C SPACING FACTOR K. THE PROGRAM COMPUTES THE AXIAL VELOCITY COMPONENTS 

C DUE TO EACH VORTEX RIRC AT A POSITION WHERE THE AIRCRAFT WING WOULD 

C BE. 4.67 DOWNSTREAM FROM THE ROTOR DISKS. THE COMPUTATION IS BEGUN 

C 33 FEET FROM THE AIRCRAFT LONGITUDINAL CENTERLINE AND PROCEEDS IN 
C I FOOT INCREMENTS UNTIL A POSITION 33 FEET ON THE OPPOSITE SIDE OF 

C THE AIRCRAFT CENTERLINE IS REACHED. THE AXIAL VELOCITY COMPONENTS 

C FROM EACH VORTEX RING BEING CALCULATED AT EAQ1 OF THE INCREMENTAL 
C STEPS. THE TOTAL VELOCITY, THE SUMMATION OF ALL THE AXIAL COMPONENTS 

C IS THEN AVERAGED BETWEEN EACH TWO SUCCESSIVE POINTS AND THIS AVERACED 

C VELOCITY IS THEN USED TO CALCULATE A FORCE (DOWNLOAD) COMPONENT. 

C THE FORMULA USED IS' F- .3*RHO*VBAft**2 . *S*C4 . THE SURFACE 

C AREA (S) IN THE FORMULA IS TAKEN TO BE A I FOOT CHORD SECTION 

C MULTIPLIED BY THE AIRCRAFT WING CHORD. 3.23 FEET. THE DOWNLOAD 

C CONFUTATION IS ONLY PERFORMED USING THE COMPUTED VELOCITIES TBAT 
C LIE BETWEEN THE ROTOR CENTERLINES. THIS>DEINC THE ONLY AREA THAT IS 
C ACTUALLY AIRCRAFT WING. THE DOWNLOAD IS THEN SUMMED AND DIVIDED BY 

C THE TOTAL ROTOR PRODUCED THRUST TO ARRIVE AT DOWNLOAD AS A PERCENT OF 

C ROTOR THRUST. THE PROGRAM IS DESIGNED TO SIMULATE A NUMBER! OF HOVER 
C CONFIGURATIONS. ANY OF THREE* THRUST COEFFICIENTS ( .60*4. .01 . .0103) 

C CAN BE USED AS ANY OF FOUR ROTOR HEIGHT ABOVE GROUND TO ROTOR DIAMETER 

C (H/D) RATIOS (H^D* INFINITY (OUT OF GROUND EFFECT. OCE). Il/D-2.0 UN GROUND 

C EFFECT. ICE). H^DU.O ICE. AND H/'D'O.S ICE). TUE PROGRAM ALSO MAKES 

C USE OF DRAG COEFFICIENTS DETERMINED FROM WIND TUNNEL TESTS IN WHICH 

C A W1NC. REPRESENTATIVE OF THE'XV-13 WING. WAS ALICNED WITH ITS CflORD 

C PERPENDICULAR TO TUE FREE STREAM VELOCITY. THE WING WAS EQUIPPED 

C WITH FLAPS. ALSO REPRESENTAT1 VE OF THE XV-13, AND TESTED SO AS TO 

C DETERMINE DRAG COEFFICIENTS CORRESPOND I NG TO DIFFERING DECREES OF FLAP 

C DEFLECTION. THE DRAC COEFFICIENTS ARE UTILIZED HERE AND HENCE GIVE 

C TOE PROGRAM THE CAPABILITY OF ANALYZING ANY OF FIVE FLAP DEFLECTION 
C CONFIGURATIONS (6.20.46.43 8 66 DECREES). 

C 

c 

Ctts*<xstsnJtss««s<»s«*«uumtsxsirs«s*»«t*igs»fst«**«xt>«««s»tunn*xnx> 

C 

C 

C 

C 

C LIST OF VARIABLES: 


C 

C CT THRUST COEFFICIENT <T/ ( RRO*A*V t**2 . > ) 

C BxD ROTOR HE! CUT ABOVE GROUND PLANE TO ROTOR 

C DIAMETER RATIO 

C FLAFSET THE DECREES OF WJNC FLAP DEFLECTION (DEC) 

C CD DRAC COEFFICIENT 

C A VORTEX RING CIRCULATION (FT**2./SEC> 

C K VORTEX RING SPACINC FACTOR 

C H ROTOR SPEED (SEOREV) 

C FID, TIME DETVEEN HOTOR BLADE PASSACES (SEC/REV) 

C KTM > CORRECTED TIME DETVEEN VORTEX niNC SHEDDING 

C <K*F( I ) ) (SEC/REV) 

C Z(l) THE POSITION OF EACH VORTEX R1NC. DOWNSTREAM OF 


a * 2 ". 
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C THE ROTOR D1ST. AS A FUNCTION OF CORRECTED TIME 

C < FT) 

C RR(I> VORTEX R IRC RADIUS (FT) (A FUNCTION OF ROTOR 

C GEOMETRY AND VORTEX RIRC DOWN STREAM POSITION 

C COORD ( I ) THE DO VR STREAM POSITION OF THE VORTEX RINCS WITH 

C RESPECT TO THE XV- IS WING POSITION < FT) 

C HOVERD THE RATIO OF ROTOR HE I CUT ABOVE CROURD PLANE 

C TO ROTOR DIAMETER (H^D) 

C NUMRINC THE HUMBER OF VORTEX RINGS. FOR AN IH GROUND 

C EFFECT (ICE) HOVEn CONDITION. THAT ARE NEEDED TO 

C 6ATISFY THE PRESCRIBED ROTOR HEICHT TO DIAMETER 

C RATIO POSITIONING OF THE CROUND PLANE. 

C Tim THE TOTAL NUMBER OF VORTEX RINCS USED IN TOE 

C ANALYSIS (DICTATED BY THE HOVER CONDITION) 

C TVICENR THE NUMBER OF THE VORTEX IllHC THAT CORRESPONDS 

C TO THE MIRROR IMAGE OF THE FIRST VORTEX RING 

C FROM THE RICITT ROTOR (ICE IJOVEn) 

C THRJEENR THE NUMBER OF THE FIRST VOnTEX RINC FROM THE LEFT 

C ROTOR (ICE HOVER) 

C RS THE TOTAL NUMBER OF VORTEX RINCS ON THE RICUT SIDE 

C (ICE HOVER) 

C LS THE NUMBER OF THE FinST VORTEX RING ON THE LEFT SIDE 

C (ICE HOVER) 

C L THE DISTANCE FROM THE AIRCRAFT LONGITUDINAL CENTER- 

C LINE TO mE POINT AT WJICI1 THE AXIAL VELOCITIES ARE 

C CALCULATED (FT) 

C RAD THE DISTANCE FROM THE AIRCRAFT LONGITUDINAL CENTER- 

C LINE TO THE ROTOR CENTERLINES (FT) 

C X THE RATIO OF VORTEX RING DOWNSTREAM POSIION TO RINC 

C RADIUS (Z(I)^RR(I)) 

C R THE RATIO OF TUE DISTANCE, FROM THE POINT AT vniCH 

C THE VELOCITIES ARE TO' BE CALCULATED TO TUE ROTOR 

C CENTERLINES. TO TUE RING'S RADII 

C S.T.U.E.B.C. CIRC. .COMPONENTS OF THE EQUATION USED TO DETERMINE THE 
C AXIAL VELOCITY COMPONENTS DUE TO EACH niNG (FROM 

C “AERODYNAMICS OF PROPULSION'. KUCIIEMANN 0 WEBER. 1953) 

C P.D INTERMEDIATE STORACE VARIABLES 

C V(M.N) THE COMPONENT OF THE AXIAL VELOCITY AT PONT H DUE TO 

C RINC N (FT/SEC) 

C VTOT(J) THE TOTAL AXIAL VELOCITY AT F.ACH OF 71 POINTS DUE TO 

C SUMMING THE COMPONENTS FROM EACH VORTEX RING (FT/SEC) 

C VBAR THE VELOCITY OBTAINED BY AVERAGING THE TOTAL AXIAL 

C VELOCITIES AT TWO SUCCESSIVE POINTS (FT/SEC) 

C DWNLD1 THE FORCE (DOWNLOAD) OBTAINED BY USING VBAR AND A 

C 5.23 SQ. FOOT WING SECTIONAL AREA (LB) 

C < DVNLDI = . 5*[UK)*VBAR**2 . *5 . 25*CI> I 

C DVNLD2 THE SUM OF TUE INDIVIDUAL DOWNLOAD COMPONENTS (LB) 

C PRCNDLI TUE PERCENTAGE OF TOTAL nOTOR THRUST THAT IS DOVN- 

C LOAD ( ( DVNLD2/TURUST )* 100 . ) < THRU ST = f (CT> ) 

C 

C 

C 

c 

I NTECER TV | CENR.TTJREENR.TNR . Y . NUHR I NC . RS . LS , FLAT SET 
REAL K.KTC72) .ZAC. LAC 

DIMENSION Z( 144). COORD < 1 44 ) . F ( 72 ) . HR ( 144) .VT0T(7I ) . V(7I . 144) 

C 

C 

C»***s*«*s****»*«*»*, ««*»,,,**»» *»,**»««»*,,,»*,,, 

C 



c 

C THE USER 19 FIRST TR Off PT LD FOR THE. DESIRED THRU ST COEFFICIENT (CT). 
C ROTOR HZ I CRT ABOVE CROUND PLANE TO ROTOR DIAMETER RATIO (H^D) . AND 

C FLAP SETTING (FLAPS ET) 

C 

C 


c*** 

c 


“****** fcUt ^" OT Ui “* 

c 

VR ITE ( 6.9) 


9 

FORMAT ( * WHAT THRUST COEFFICIENT (CT) WILL 
WRITE ( 6 . 1 1 ) 

BE DSED7* ) 

1 1 

FORflATC * ( .0094, .01 . .0105)7* ) 
READ(5, 10)CT 


10 

FORMAT ( C 1 2 . 3 ) 
WRITE( 6 *16) 


1 6 

FORMAT ( * WHAT H/D RATIO WILL BE USED? ’ ) 
VR!TE(6, 14) 


14 

FORMAT ( ’ (0.5. 1 .0,2;0. OR OCE HOVERD’ ) 
VRITE<6, 17) 


17 

FORMAT(* (ENTER 0.0 FOR OCE HOVER)* ) 
READ ( 5 . 1 81H0VERD 


IB 

F0RMAT(C12.3) 
VRITE(6, 19) 


19 

FORMAT ( * WHAT FLAP SCTTINC WILL DE USED?* ) 
WRITE! 6.13) 


18 

FORMATC (0,20.40.43,60 DEC.)*) 
READ (5. 8) FLAPSET 


Q 

FORMAT (12) 


C 



C 



c*«* 

c 

c 

c 

c 

c 

FROM THE SELECTED FLAP SETTING. THE PROPER 

uttttstxuvt'vssitt'us'Jisutxaji 

DRAC COEFFICIENT IS SET 


c 

c 

1 F ( FLAPSET . EQ . 0 > THEN 

CD- I .70 

ENDIF 

IF(FLAPSET . LQL 20 1 THEN 

CD* 1 .52 

ENDIF 

I F( FLA PS ET.BQ. 40) THEN 

CD* I .32 

ENDIF 

IF (FLAPSET. EG, 45) THEN 

CD* 1 .26 

ENDIF 

I F( FLAPSET. EO. 60) THEN 
CD 3 I .08 
ENDIF 
C 
C 

OixmimtitimititiimtituitmnuttstttfratiittimsnrxfitnAtxntffn 

c 

c 

C THE VALUES FOR VORTEX RINC CIRCULATION (A) AND VORTEX RING SPAC1NC 



92 


c factor <k> arc how srr corresponding toithe selected t hhu st i coot i c i err 



c 

IF<CT.EQ. .0494) THEN 

4009.0a 
1C* *446 
EHDIF 

IFCCT. EQ. .01 ) THEN 
A* 20 1 . 1 6 
K* .446 
EHDIF 

IFICT.EQ. .0103) TBEH 
A-211 .22 
K* .442 
EHDIF 

C 

c 

C*»*** ****** ******** 

C 

C 

C SINCE THE AR I CRAFT HAS A THREE B UVD ED ROTOR, TTUIEE VORTEX RIHCS ARE 

C CONT I RUOUSLV SHED FROM THE ROTOR IN A HELICAL SHAPE. AN APPROXIMATION 

C IS HADE HERE THAT TRANSFORMS THE THREE CONTINUOUS HELICAL VORTICES 

C INTO THREE ROUND VORTEX RINGS. EVENLY SPACED IN TIKE AND SPACE. OVd 
C ONE ROTOR REVOLUTION. THE VORTEX RINGS’ POSITION IN TIKE IS THEN 

C ALTERED BY THE APPROPRIATE VORTEX RING SPACING FACTOR <K> THAT WAS 

C PREDETERMINED USINClTBE CALIBRATION PROGRAM VERSION 

C 
C 

C*mtm«t**m*i>xnm»»><«tt>*miitx«ntitn>n«txu>ttiatxxxx*<m>xn 

C 

C 

H* . 1022 
DO 20 1*1,72 

F < I )* l*B^3.0 
KT( I ) *K*F< I 1 
26 CONTINUE 

C 
C 

c************* ******************************** ********************* ***** ** **~ * 
c 

C 

C DEPENDING ON WHICH THRUST COEFFICIENT WAS SELECTED. CONTROL IS NOW 
C TRANSFERRED TO AN APPROPRIATE FIFTH ORDER INTERPOLATING POLYNOMIAL. 

C THE POLYNOMIALS ARE. USED TO SET THE DOWNSTREAM POSITION OF EACH VORTEX 

C RING AS A FUNCTION OF THEIR RESPECTIVE CORRECTED TIKE SPACING. THE 
C RADII OF EACH VORTEX RING IS ALSO CALCULATED. THESE BEINC A' FUNCTION 

C OF THE RING’S DOWNSTREAM POSITION AS WELL AS THE ROTOR’S DIAMETER. 

C THE POSITION OF EACH RINC WITH RESPECT TO THE AIRCRAFT WING* IS ALSO 
C DETERMINED. 

C IF ONE OF THE ICE BOTOV CONDITIONS 19 DEINC ANALYZED. CONTROL IS 

C TRANSFERRED TO ANOTHER SECTION OF THE PIVOCRAM WllUf THE DOWNSTREAM 

C POSITION OF A VORTEX RINC CLOSELY. IF NOT EXACTLY. SATISFIES THE 

C PRESCRIBED H/D RATIO. 

C 

c 

CxmixixxixxixixxxiixiiixxiixmiiiiiixxiimtiitiiixmxxxmiixiitixtmtM 

C 
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C 

IF(CT.EQ. .01 ) CO TO 30 
(F(CT.£a. .6103) CO TO 33 
C 
C 

23 DO 40 1*1.36 

Z< I >*“63 .835*KT < 1 >**5.4208 . 424*KT( 1 )**4 . -242 . 599*KT ( I )**3. 

♦ ♦120.79O»KT ( 1 >**2. +73. 9)2*KT( I > 

RR( I >*0.8373*SORT( 1 .✓( . 5* < 2 . -EXP (- .643 1 2*Z ( !>>)>) 

COORD ( I )*Z( I >-4.67 

1F( (ROVERD.EO.e.S) .AND. (Z( I ) . CE . 12 . 5 > . OR . 

♦ (HOVERD. EQ. 1 .0) .AND. (Z( I > .CE , 23, > .OR . 

♦ (ROVERD.EQ.2.0) .AND. <Z( l > .CE.50. > > THEN 
NUttRI NC= I 

ZAC=Z( 1 ) 

CO TO 26 
END1F 
46 CONTINUE 

CO TO 24 
C 
C 
C 

30 DO 43 1=1 .36 

Z ( I )s-60. I52*KT( I ) **5.4221 .77*KT( l >**4 . -236 .836*KT< I )**3. 

♦ +126.62*KT( I )**2.478.0G3*KT< I ) 

RR( 1 > *8.837 5* SORT ( l ./( .5*12. -EXPt “ . 09512*Z( !>))>> 

COORD ( I ) «Z( 1 >—4l 67 

I F ( < HOVERD . EQ .0.3) . AND . < Z< 1 > .CE. 12.5) .OR. 

4 (HOVERD. EQ. 1 .0). AND. <Z< 1 > .CE .23 . ) .OR . 

4 (HOVERD. £0.2.0) .AND. (Z( l ) . CE. 50 . > ) THEN 
NUT1R1NC*! 

ZAC*Z< I ) 

CO TO 26 
ENDIF 

43 CONTINUE 

CO TO 24 
C 
C 
C 

33 DO 50 Ial.36 

Z( I ) * “64 .7 46*KT ( I >**3.421 | .756*KT< I > **4 . “247 . 2 1 9*KT < I >**3 . 

♦ >I23.687*KT( I >**2.40! ,398*KT( l ) 

RR( I >=8.8375*SQRT( 1 ./( . 5* (2 . -EXP ( - . 093 1 2*Z< !)>)>) 

COOWH I > *Z( I >-4.67 

I F < (HOVERD. EO. 0.3) .AND. <Z( I > .CE. 12 . 3) .OR. 

♦ ( HOVERD. EO . 1 .0) .AND. (Z(l > .CE.25 . > .OR . 

4 (HOVERD. Ett. 2.0) .AND. (Z( 1 ) .CE. 30 .) ) THEN 

NimRINC* I 
ZAC=Z( I ) 

CO TO 26 
ENDIF 

30 CONTINUE 

C 

C 

Cmxvn<msntti*num«*txs»s>tSMxikt»«stxsi*t>*si>**nusstnun«*** 

C 

C 

C IF AN OCE BOVER CONDITION IS DEINC ANALYZED. THIS SECTION PRODUCES 

C THE VORTEX RINC SYSTEM FOR THE LEFT ROTOR. THIS LF.FT UOTOR VORTEX 

C RINC SYSTEM 19 IDENTICAL TO THE RICHT EXCEPT FOR TUE FACT THAT ITS 

C CENTERLINE COINCIDES WITH TEAT OF TUE LEFT ROTun . CONTROL IS TUEN 



ORIGINAL PAGE IS 
OF POOR QUALITY 


C TRANSFERRED TO STATTHERT LABEL 60. 

C 

C 

c 

c 

24 DO 41 1>I .36 

Zt36+I >*Z< I > 

RJU36+I )=RR( I ) 

COORD ( 36+ 1 ) x COORD ( I ) 

41 CONTINUE 

TNR=72 
CO TO 66 
C 
C 

C»>tmf<iimni*x«t***»xmt«xmtt«m«ttxtuyu«snmxss«*nxmmn>n 

C 

c 

C IF AN ICC nOVCR CONDITION IS BEING ANALYZED. CONTROL IS TRANSFERRED 

C TO THIS SECTION FROM THE INTERPOLATING POLYNOMIALS. THIS SECTION 

C PRODUCES THE MIRROR* IMAGE OF THE VORTEX RINC SYSTEM ON THE RICIfT 

C SIDE. REFLECTED ABOUT THE GROUND PLANE DICTATED BY TOE fl/D RATIO. 

C ONCE THIS IS ACCOMPLISHED. AN) IDENTICAL VORTEX RING SYSTEM IS PRO- 
C DUCED FOR TOE LEFT SIDE. 

C THE FIRST SECTION SUPERIMPOSES A VORTEX RING. IDENTICAL IN RADIUS 
C AND POSITION. ONTO THE LAST RING THAT WAS PRESCRIBED FROM TOE INTER- 
C POLATINC POLYNOMIAL. THE RING THAT SATISFIES THE B/D RATIO. 

C 

C 

Cxxuxxxx*xxx*kxx*xxxxx**xx*x*s*s«******x*x*x**x*xx*****x* ******************** 

C 

C 

26 Y- NUMRINC+ I 
Z(Y)-Z(NUMRINC) 

RR( Y ) *RR< NUMR1 NC) 

COORD < Y )* COORD < NUMR I NC ) 

C 

C 

C *********************************** *** *»* **■*•« ****T***^««** *** ***»*'«y«*** : »*»** 

c 

c 

C THIS DO LOOP PRODUCES THE REMAINING * MIRRORED* VORTEX RINC SYSTEM- - 
C BECINNINC AFTER THE TWO SUPERIMPOSED VORTEX RINC S. THE FIRST RINC 
C PRODUCED BY THIS SECTION IS IDENTICAL IN RADIUS AND POSITION. WITH 

C RESPECT TO THE SUPERIMPOSED R1NCS, AS TOE RINC THAT OCCTJRRS PRIOR 

C TO THE SUPERIMPOSED RINCS. HENCE. THIS SECTION PRODUCES A SET OF 

C VORTEX RINCS THAT ARE IDENTICAL IN RADIUS AND POSITION WITH ; RESPECT 

C TO THE SUPERIMPOSED RINCS. 1 . o . • CnOUND PLANE. 

C 

C 

C*»*nxx*s»xxxxxxXxxxxxsxa*xxxxxxx*xxx*»xxxXxr**x**xYxx*x»x=**x*xxxxxxxxsxxxx* 

C 

c 

DO 27 I *2. NUMRI NC 
LOPE* I ♦NUMRI NC 
LI P* Y-l 

Z< LOPE) *21 Y) ♦IZIY ) -ZIL1P ) ) 

RR< LOPE) *RR( LIP) 

COORD ( LOPE ) * Z ( LOPE ) -4 . 67 

27 CONTINUE 
C 
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C 

C 

C 

C THIS SECTION PRODUCES A VORTEX RING SYSTEM FOR THE LEFT SIDE THAT 
C IS IDENTICAL TO THAT PREVIOUSLY PRODUCED FOR THE RIGHT SIDE- 

C 
C 

C 

C 

TW1CENR«2*NUMR1NC 
THREENR*3*NUMR I NC 
TNR=4*NUHRINC 
DO 28 I»1 .TVICENR 

ZM*TWICENR)=ZM > 

RAM ♦TVICENR) S RR< I ) 

COORD C I ♦TV I CEMR ) » COORD ( 1 ) 

28 CONTINUE 

C 

C 

C»m«n«as»iukzt>xut»i**nsist«ixst»»»mx»x»xxxstxutx«sxmm» 

C 

C 

C THIS SECTION SETS UP A CUTOFF 1 POINT REFERRING TO THE NUMBER OF 

C VORTEX RINCS THAT OCCUR ON THE RICIIT AND LEFT SIDE OF THE SYSTEM. 

C THESE CUTOFF POINTS ARE LATER USED TO DETERMINE THE DISTANCES 
C BETVEEN POINTS AT WHICH THE AXIAL VELOCITIES ARE CALCULATED AND 
C THE APPROPRIATE VORTEX SYSTEM CENTERLINE. 

C 

C 

Cstm*m»»*ci*<»xsxxt»ss*»sxi<>***ssumxtixussss»xt«s»amsmm>** 

C 

C 

60 IFMNR.EQ.36) THEN 

RS= 18 
LS= 19 
ENDIF 

IF<TNR.EQ.4e> THEN 
RS = 2e 
LS*2 I 
ENDIF 

IF<TNR.E0.68> THEN 

RS*34 

LS=33 

ENDIF 

1F(TNR.EO-72) THEN 

RS*36 

LS*37 

ENDIF 

IF<TNR.EQ. 128) THEN 

RS*64 

LS»65 

ENDIF 

I F < TNR . £0 . 1 32 ) THEN- 

RS*66 

LS= 67 

ENDIF 

IFlTNR.EO. 136) THEN 

RS=68 

LS* 69 
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ENDIF* 



THIS SECT I ON IS THAT IN VHIOTlTHE ACTUAL AXIAL VELOCITY COMPONENTS 
ARE CALCULATED. THE SECTION CONSISTS OF A SERIES. OF 1F/THEN STATE- 
MENTS AND TWO NESTED DO LOOPS'. 

THE PURPOSE OF THE IF/THEN STATEMENTS IS TO SET THE CORRECT RELATION- 
SHIPS REGARD INC THE' DISTANCES THAT ARE NEEDED FOR: THE VELOCITY CAL- 
CULATIONS. 

THE CALCULATION BEGINS AT A POINT 35 FEET FROM THE AIRCRAFT LONGITUDINAL: 
CENTERLINE AND ENDS AT A POINT 33 FEET ON THE OPPOSITE SIDE OF THE 
CENTERLINE. THE IF/THEN STATEMENTS SET THE PROPER DISTANCES BETWEEN 
EACH SUCCESSIVE POINT. WHERE THE VELOCITIES ARE TO BE CALCULATED, AND 
THE CENTERLINES OF THE TWO RESPECTIVE VORTEX SYSTEM CENTERLINES 
(♦/-l 6. 6833 FT). THE DISTANCES ARE ALSO DEPENDENT ON WHICH SET OF 
VORTEX RINGS ARE BEING CONSIDERED RELATIVE TO THE LOCATION OF THE 
POINT AT WHICH THE VELOCITIES ARE BEING CALCULATED. 

WITH THE PROPER RELATIVE DISTANCES ESTABLISHED. THE RATIOS X <RINC 
DOWNSTREAM POSITION TO RING RADIUS) AND R I DISTANCE. FROM THE PO INT TO 
THE CURRENT VORTEX RING SYSTEM CENTERLINE. TO RINC RADIUS) ARE DETER- 
MINED AND THE VELOCITY CALCULATION IS CARRIED OUT. 

IF AN ICE HOVER CONDITION IS BEING ANALYZED. AN IF/THEN STATEMENT 
CHANCES THE CIRCULATION SO AS TO MAKE IT NEGATIVE IN THE REFLECTED 
VORTEX SYSTEM. 

FINAL VELOCITY RESULTS ARE STORED IN THE DOUBLEY SUBSCRIPTED ARRAY 
V(H,N) , ACCOMODATINC VELOCITIES FROM 71 POINTS AND TNR VORTEX RINGS. 




L«-33 

DO 63 H* I .71 

DO 76 NM .TNR 
RAD 8 t 6 . 6833 

IFC (L.LT.0) .AND. (N.LE.RS) ) THEN 
J* -L 

RAD 8 16.6833 
END I F 

IFUL.LE.-m.AND. IN.CE.LS> ) THEN 
J*-L 

RAD 8 -I 6.0833 
ENDIF 

I F ( ( L. CE .-16). AND. (L.LT.O) -AND- (N.CE.LS) > THEN 
J = L 

RAD* 16.6833 
ENDIF 

1FML.CT.6) .AND.CN.CE.LS) > THEN 
J = L 

RAD* 1 6.6833 
ENDIF 

IFC < L. LE. 16) .AND. ( L. CT . 6 ) . AND. ( N . LE. RS ) ) THEN 
J*-L 

RAD- I 6.9833 

ENDIF 

IFC (L.CE. 17) .AND. ( N.LE.RS) ) THEN 
J*L 
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RAD«- 1 6.0833 
END1F 

(F<L.Ca.O) T HEN 
J*L 

rUO= 16.0833 
END IF 

X*-COORD<N)'RR<N> 

R* (RAO+J >/RR(N> 

S“ < 4 . *R>^<X**2. -MR* I . )**2 . ) 

T*ALOC ( 4 . ✓ < I . -S >** . 5 ) 

U=T+< (T-l . )/ 4 . >*( I .-S>*< 9./f»4. >*(T-<7./6. ) )*< 1 .-S)**2. 

EM . ■*■( . 3* (T- . 5 l * ( 1 . ) > ♦ ( 3 . ✓ I 6 . I*< T— < 1 3 . / 1 2 . ) ) * < I . -S>**2. 

B= 1 ./IX«2. + (R*I . )**2. >**.5 

C=-< I .♦( <2.*<R-I . ) )^<X**2. + <R-I . )**2 . ) > > 

D = tl+C*E 

P S B*D 

A-ADS(A) 

I F < (UOVERD.NE.O.0) .AND. ( N .CT. NUMRl NC > .AND. ( N . LE .TV I CENR ) 

♦ .on. 

♦ (I10VEJID . NE.O - 0) . AND . < N . CT.THREENR ) . AND . ( N . LE . TNR ) > THEN 
A»-A 

END IF 

CinC = A/(2.*3 . l4l39*fUUN> ) 

v<n.N )=cinc*p 
70 CONTINUE 

C 

L 5 L+ 1 

63 CONTINUE 




THE NEXT NESTED DO LOOPS PERFORM THE TASK OF SUTOIINC T1IE VELOCITIES 
AT EACH POINT DUE TO EACH OF THE !l I NOS . IIKNCC. FROM I1LHE. THE TOTAL 
VELOCITIES AT EACH OF 7 1 POINTS ARE OBTAINED AND STORED IN THE AJUUY 
VTOT ( J ) . 


timxisxmuyiimxftttmxnifSixittxmytutmmtmmmuiLniii ******* 


DO 75 J- 1 ,71 

DO 76 1 = 1 ,TN7l 

VTOT < J ) * VTOT ( J ) ♦ V ( J . I ) 

6 CONTINUE 

VR1TEI * , 88 > J . VTOT ( J ) 

8 FORMAT ( s . 2X . 5HVT0T ( . J2.2H ) = . F6 .2 .3X . 60FT/SEC ) 

5 CONTINUE 


THIS SECTION OUTPUTS THE CITOSEN THRUST COEFFICIENT (CT). ROTOn nEICHT 
ABOVE GROUND TO KOTOR DIAMETER RATIO (O/Dt. THE ASSIGNED VALUE OF THE 
VORTEX RING SPAC1NC FACTOR <K) AND THE RESULTING NUMBER OF VORTEX RINCS 
<TNR) . 
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c 

c 

LAC* 2 3 . 

HACTRD* ZAO’ LAC 

WR I T E ( 6 . 3 3 ) CT . HA CTRD 

33 FORttATt//.* FOR CT* * .F6.5.* ,A TRUE H/D«* .F4.2) 

I F < BOVERD . EQ .0.0) THE* 

WRITE! 6 . 34 ) 

34 FORMATC <B/D«0.00 SIGNIFIES AN OCE HOVER CONDITION > * 1 
END IF 

WniTE( 6 . 83 )IC ,TNR 

83 FORMAT! ’ K= * . F4 . 3 . * AND*. IX. 13.* VORTEX R t NCS * ) 

C 

C 

c 

c 

C THE FOLLOWING I>0 LOOP PERFORMS A SERIES OF TASKS. FIRST, ONLY TT1E 
C COMPUTED VELOCITIES THAT LIE BETWEEN THE ROTOR CENTERLINES ARE CONDIDERED. 

C THIS BEINC WHERE THE ACTUAL AIRCRAFT WING IS. START INC ON THE LEFT 

C HAND SIDE. EACH TWO SUCCESSIVE VELOCITIES ARE AVERAGED. TUE AVERACED 

C VELOCITIES THEN DE1NC USED TO COMPUTE A DOWNLOAD. IF THE AREAS WHERE 
C THE VELOCITIES OCCUR ARE IN THE VICINITY OF THE WINCS AND NOT THE 

C FUSELACE. THE PREDETERMINED DRAG COEFFICIENT IS USED. IF THE AREAS 

C WHERE THE VELOCITIES occun ARE IN THE VICINITY OF THE FUSELAGE. THE 

C DRAC COEFFICIENT IS ASSUMED TO BE I. FOR WANT OF A BETTER VALUE. 

C TWO 1F/TOEN STATEMENTS TAKE CARE OF THIS DISTINCTION. THE DOWNLOAD 

C COMPONENTS ARE THEN SUMMED TO ARRIVE AT A TOTAL DOWNLOAD. THIS FINAL 

C VALUE (IN POUNDS) IS THEN OUTPUT ALONC WITH TUE CHOSEN FLAP SETTING. 

C 

C 

C 

C 

DO 699 1-20.52 

VBAR= ( VTOT( I ) ♦VTOT ( t + l > )/2.0 
I F < ( I .CE.33) .AND. ( I -LE.39) ) THEN 
DVNLDl e . 5* . 002370* VBAR* *2 . 0*5 . 23 
ENDIF 

1 F ( ( I . LT . 33 ) . OR . < I . CT . 39 ) ) THEN 
DVNLDl <= .5* .002378»VBAi\**2.0*5.23*CD 
END l E 

D WN LD2 * D WN LD2 ♦ DWN U) 1 
6>99 CONTINUE 

VR1TE( 6 .399 1DVNLD2 

399 FORMAT <* THE XV- l 3 ROTOR DOWNLOAD IS «’,F0.3.* LD . * ) 

VR | TE( 6 . 589 ) FLAPSET .CD 

389 FORMAT ( ' FORA FLAP SETTI NC OF * . I X . 1 3 . ’ DECREES ( Cd - * . F4 . 2 , * ) . * ) 

C 

c 

C 

C 

C TRE LAST SECTION CONSISTS OF TTTREF. I F'TTTF.N STATEMENTS AND RANDLES 

C THE TASK OF CALCULATING AND OUTPUTTING THE FINAL TOTAL DOWNLOAD AS 

C A PERCENTAGE OF TOTAL ROTOR THRUST. SINCE THE TOTAL ROTOR THRUST 

C PRODUCED IS A FUNCTION OF THE 1 THRUST COEFFICI ENT . TRE PREREQUISITE 

C FOR THE IF/THEN STATEMENTS IS THE THRUST COEFFICIENT. EACH: IF-THEN 

C STATEMENT DIVIDES THE TOTAL DOWNLOAD BY TUE APPROPRIATE TOTAL THRUST 

C AND OUTPUTS THIS PERCENTAGE VALUE. 

C 



99 

ORIGINAL PAGE IS 
OF POOR QUALITY 




C 

1FCCT.EQ. .3*94) THEN 

PRCKDLl * < DWNL02/ 12956.621*1 00 .• 

WRITE! 6 , 400 > 

400 FORMAT! * BASED ON A THRUST COEFFICIENT OF .0094’ l 
VRITE(6.40I JPRCKDLI 

401 FORMAT!* TUIS REPRESENTS’ .IX.F6.3,* r. OF TOTAL nOTOR THRUST 1 . ) 
EHDIF 

C 

IFICT.EQ. .01 > THEN 

PRCNDLI = ! DWNLD2/ I 3703 .641*100. 

VRITE( 6 . 403 ) 

403 FORMAT! ’ BASED OR A THRUST COEFFICIEITT OF .01' 1 
VR|TE( 6.404 ) PRCNDL 1 

404 FORMAT ( ’ THIS REPRESENTS * , I X . F6 . 3 . * n OF TOTAL ROTOR THRUST’ > 
END IF 

C 

IFfCT.EQ. .0103) TIIEU 
PRCNDLI MDVNLD2/1 4472. 02 1*1 00. 

WRITE! 6 ,406 ) 

406 FORMAT! ’ BASED OR A THRUST COEFFICIENT OF .0105’ ) 

WRITE! 6 . 407 1 PRCNDL l 

407 FORMAT!* THIS REPRESENTS* , 1 X . F 6 . 3 . * 7. OF TOTAL ROTOR THRUST 1 . ) 
ENDIF 

C 

STOP 

EMD 



Appendix F 

Calibration Version of Model Rotor 
Computer Program 
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THIS PROGRAM IS DESIGNED TO SIMULATE THE MODEL ROTOR SYST E M TEST 
PERFORMANCE THROUCH THE USE OF A VORTEX RING ANALYSIS. THE PROGRAM 
IS BASED ON A ROTOR THRUST COEFFICIENT OF .0073 AND TIP SPEED OF 
83.23 FT/SEC. THIS PARTICULAR PROCRAM' VERSION WAS DES1CNED TO : 
CALIBRATE THE VORTEX RINC SPACINC FACTOR (K>. THE VALUE FOR K IS 
AN INPUT PARAMETER IN THIS PROGRAM AND- 18 VARIED UNTIL AN AVERAGE 
INFLOW VELOCITY THROUCH THE ROTOR DISK THAT AGREES WITH THAT 
PREDICTED BY MOMENTUM THEORY i IS OBTAINED. THE POINTS AT WHICH THE 
VELOCITIES ARE CALCULATED LIE ALONG A LINE WHICH; IS COINCIDENTAL 
WITH THE ROTOR PLANE. 




LIST OF VARIABLES*.* 


K VORTEX RINC SPACINC FACTOR 

H ROTOR SPEED (SEC/REV) 

F( I ) TIME BETWEEN VORTEX RINGS (SEC) 

ETUI CORRECTED TIME BETWEEN VORTEX RINCS <SECJ 

CK*F<I>> 

Z II) THE POSITION OF EACH VORTEX RINC. DOWN- 

STREAM OF THE ROTOR DISK, AS A FUNCTION OF 
CORRECTED TIME (FT) 

RR(I) VORTEX RINC RADIUS (FT) 

(A FUNCTION OF ROTOR CEOMETRY AND RINC 
DOWNSTREAM’ POS 1 T I ON ) 

TNR THE TOTAL NUMBER OF VORTEX RINCS CONSIDERED 

IN THE ANALYSIS 

A THE VORTEX) R I NC CIRCULATION (FT**2/SEC) • 

CT THRUST COEFFICIENT (THRUST/ (RMO* AREA* V i | p**2 ) > 

POINT THE DISTANCE BETWEEN THE ROTOR AXIS AND THE 

POINT WHERE THE AXIAL VELOCITIES ARE BEING 
CALCULATED (FT) 

RAD THE ABSOLUTE VALUE OF POINT (FT) 

X THE RATIO OF VORTEX RING DOWNSTREAM 

POSITION TO RINC RADIUS <-Z<I>/RRU)> 

R THE RATIO OF POINT TO THE RINC RADIUS 

<RAD/RR( I )> 


8, T,U,E.B,C, CIRC. .COMPONENTS OF THE EQUATION USED TO 


DETERMINE THE AXIAL VELOCITY DUE TO EACH 
RING < FROMJ AERODYNAMICS OF PROPULSION* 
EUCHEMANN S WEBER. 1953) 

P.D INTERMEDIATE STORAGE) VARIABLES 

V(M.N) THE COMPONENT OF THE AXIAL VELOCITY AT 

POINT M DUE TO RINC N (FT/SEC) 

VTOT ( I ) . . THE TOTAL AXIAL VELOCITY AT EACH OF 25 

POINTS DUE! TO SUMHINC THE COMPONENTS FROM 
EACH VORTEX RINC (FT/ SEC) 

VAVG THE SUM OFiTHE AXIAL VELOCITIES ALONC THE 

ROTOR DISK (FT/SEC) 

VDIST THE AVERACE INFLOW VELOCITY ALONC THE ROTOR 

DISK I FT/SEC) 
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c 

c 

C 

C 

REAL K.KT<72) 

INTECER TTfR 

DIMENSION Z(144) ,F<36) .RR< 144) .VT0T<71 ) .V(71 . 144) 

C 

C 

CJraMaaaaaaMaaaaMaaaaaaaaxrMmaw**************-**************************** 

C 

c 

C THE VALUE TO BE USED FOR THEiVORTEX RINC SPACINC FACTOR (K) IS 

C FIRST READ IN 

C 

C 

C 

VR!TE(6.4) 

4 FORMATS VHAT K VALUE WILL BE USED?* ) 

READ ( 3 , 3 )K 
3 FORMAT <F8 . 7 ) 

C 

C 

C 

C 

C NEXT. BASED ON A ROTOR SPEED OF 700 RPK. THE TIKE SPACINC OF 

C THE VORTEX RINCS IS CALCULATED USINC THE RING SPACINC FACTOR <K) . 

C 

C 

C<y*xs<ss**itxit»iixfy>»*ss**«st*«sx*«*<xiixitx««*>sx*>«*stx*y**ss<ssx*u»xss< 

C 

C 

H = .0037 
DO 20 1-1.72 

F < I) * l*H/2.0 
KT( I )=K*F< I ) 

20 CONTINUE 

C 

C 

Cxxxxxxxrxxxxxxxx>»xxx»»txx»xxxxx»«i»j^wtxx»f»xT»xx«xxxxyxxx»x»»»rxqx»r>»> 

C 

C 

C A FIFTH ORDER INTERPOLATING POLYNOMIAL IS USED TO DETERMINE THE 

C DOWNSTREAM POSITION OF EACH VORTEX RING AS A FUNCTION OF TIME. 

C EACH VORTEX RING'S RADIUS 18 ALSO CALCULATED AS A FUNCTION OF 

C RING DOWNSTREAM POSITION. 

C 

c 

Cmxxxmxxxxxuxxxxxxtxxsxxxxxxixxximmxminnuxxxxxitximxnuxxst 

C 

c 

23 DO 40 1*1 .72 

Z<I ) *3 , 342 1 »KT( I )a*3. -M 3-. 3269*JCT( J ) **4 . -24 . 81 07*KT< 1 ) 

♦ **3 . +12 .B299*KT< I ) **2 . *7 . 2430*KT ( I ) 

RR< 1 )» .80274*SORT< 1 .✓( . 3a <2 . -EXP ( - I .0476*21 !)>>)) 

40 CONTINUE 

TNR*72 . 

C 



103 


ORIGINAL PAGE IS 
OF POOR QUALITY 


THE VALUES OF THRUST COEfTICJEHT <CT> . RISC CIRCULATION <A) AJTO THE 1 
POSITION OP THE FIRST POIRT. WITH RESPECT TO THE) ROTOR AXIS. (POINT):. 
ARE NOV ENTERED. 


m wmmm *tt *9nm*.** M*wm** *** 


CT= .0073 
A-2.226 
POINT* -1 .3 


TWO DO LOOPS ARE NOV USED TO CALCULATE EACH VORTEX RING’S AXIAL 
VELOCITY COWPOirrarr AT each point of interest, the inner loop 
CALCULATES THE VELOCITIES FROM EACH OF THE RINGS AND THE OUTER 
LOOP INCREMENTS THE POINT AT WHICH THE VELOCITIES ARE CALCULATED. 


9m***M***x*n-*****m*** n***** n*M*xm*M****n ** m * *Mm* 


DO 63 H*1 .23 


DO 70 N* f ,T7CR 

I Ft POI NT . LT .0 . ) RAD«*POINT 

iFtPonrr.CT.e. ) rad*point 
X*-Z(N)/RR(N) 

R-RAD/RR(N) 

S*(4.*RI/(X«2.-MRM . >»* 2 . ) 

T * ALOC ( 4 . /SORT I 1 . -S> > 

U = T+< IT- 1 . )X4. >*U .-S)*< 9.^64. )*CT-<7./6. ) )*< I . -9>*»2. 

EM .♦( . 3* <T- . 5 )* < l . -S))-M3./T6. )* (T-( 13./|2.))*<1 .-S>«2. 
B- 1 . /SORT ( X**2 . ♦ ( R* I : )**2 . ) 

C*-( I .♦< (2.*(R-| . ) )/<X**2 . MR-1 . >**2. ) ) ) 

D»U4C#E 
P* B*D 

CIRC- A/ (2. *3 . 14t56*RR(!U ) 

V<H.N>*CIRC*r 

CONTINUE 


POINT-POINT*. 123 
CONTINUE 




THE VELOCITY COKPONTNTS DUE TO EACH OP' THE VORTEX RINGS ARE NOV 
SUMMED TO Ail I VE AT! A TOTAL AXIAL VELOCITY AT EACH OF THE 23 POINTS 
OF INTEREST. 
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c 

C**m ************ ************ &** » m>imm»imi«*i ttww nmnnummi m 

C 

C 

DO 73 J«1 .23 
C 
C 

DO 60 1 = 1 .TNR 

VT0T ( J ) * VTOT ( J ) *V < J . 1 ) 

ee continue 

c 

c 

73 continue 

c 

c 

C 

C 

C THE TOTAL AXIAL VELOCITY AT EACH POINT OF INTEREST IS NOV PRINTED 

C OUT ALONG WITH THE ASSOCIATED VALUES OF THRUST COEFFICIENT <CT). 

C VORTEX RINC SPACINC FACTOR AND THE TOTAL NUMBER OF VORTEX RINGS 

C USED IN THE VELOCITY CALCULATIONS . 

C 

C 

C 

C 

VRITEC 6,33 >CT 

33 FORMAT! * FOR A THRUST OEFFICIENT OFCT=* ,F6.5) 

VRITEC6 ,83 JK.TNR 

83 FORMAT! IX , * K«*.F8.7.* AND*, IX. 12.* VORTEX R1NCS' ) 

VRITE( 6 . 90 > 

FORMAT! IX.* THE TOTAL AXIAL VELOCITY AT EACH POINT ALONC A» LINE* ) 

WRITE! 6,93) 

93 FORMATUX,’ COINCIDENTAL WITH TUE PLANE OF THE ROTOR IS:*) 

C 

C 

DO 97 ,25 

VR ITE < 6.96) J.VTOT(J) 

96 FORMAT! /.2X,3HVT0T< .12,23)* . Cl 2 . 5 . 6HFT/SEC ) 

97 CONTINUE 
C 

c 

C 

c 

C SINCE FOR CALIBRATION PURPOSES THE AVERACE VELOCITY ALONC THE 

C ROTOR DISK IS WHATS OF PRIME IMPORTANCE. THAT VALUE IS NOV 

C CALCULATED AND OUTPUT. 

C 

c 

Cn*tmatuuwujtstMt<stmms*iM>t>»*i>s*xitt«mia«x»aMi*mxxsntua 

C 

c 

do tee t=4.i3 

VAVC* VAVC+VTOT ( I ) 
iee CONTINUE 

c 

VDiSK*VAvc/ie .e 

WRITE ( 6 , 1 16 JVDISK 

lie FORMAT (x , * THE AVERACE JNFLOV VELOCITY IS* , 1 X. FO. 5 , * FT /SEC*) 
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C 

c 

C THIS PROGRAM IS DESICNED TO SIMULATE THE MODEL ROTOR SYSTEM TEST 

C PERFORMANCE THROUGH THE USE OF A VORTEX RINC ANALYSIS. THE PROGRAM 

C IS BASED OR A ROTOR THRUST COEFF I C I ENT OF .0073 ARD TIP SPEED OF 

C 63.23 fT/SEC. UTILIZING THE1 PREDETERMINED VORTEX RIRC BPACIRC 

C FACTOR K. THIS PROGRAM COMPUTES A DO VR LOAD BASED* OR THE MODEL VI RC 

C (CHORD ■ 4.7 IR. . SPAR « 16.81 IR.) AT A POSITIOR 1.5 IN. DOWN- 

C STREAM OF THE ROTOR PLANE. THE COMPUTED AXIAL VELOCITIES ARE USED 

C TO DETERMINE A VALUE OF DYNAMIC PRESSURE <q>. VBICB WHEN MULTIPLIED 

C BY A SECTIONAL SURFACE AREA CIVES A DOWNLOAD COMPONENT. 

C 

c 


C 

C 




C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


LIST OF VARIABLES: 


K VORTEX RING SPACING FACTOR 

H ROTOR SPEED (SEC/REV) 

F( I ) TIME BETWEEN nOTOR BLADE PASSACES (SEC/REV) 

KT < I ) CORRECTED TIME BETWEEN VORTEX RINC SHEDDING 

<K*F ( I ) ) (SEC/REV) 

Z(I) THE POSITION OF EACH VORTEX RINC. DOWN- 

STREAM OF THE ROTOR DISK. AS A FUNCTION OF 
CORRECTED TIME (FT) 

COORD ( I ) THE DOWNSTREAM POSITION OF TT1E J VORTEX RINCS 

WITH RESPECT TO A LINE THAT IS COINCIDENTAL 
WITH THE MODEL ROTOR’S VI NC (FT) 

RR(1) VORTEX RINC RADIUS (FT) 

(A FUNCTION OF ROTOR CEOMETRY AND VORTEX RINC 
DOWNSTREAM’ POSITION) 

NUMRINC THE NUMBER) OF VORTEX RINCS. FOR AN ICE ROVER 

CONDITION. THAT SATISFIES THE PRESCRIBED ROTOR BE1CHT 
ABOVE GROUND PLANE TO ROTOR DIAMETER RATIO 

TKR THE TOTAL NUMBER OF VORTEX RINCS CONSIDERED 

IN THE ANALYSIS 

A THE VORTEX) RING CIRCULATION <FT**2/SEC) 

CT THRUST COEFFICIENT (THRUST/ ( RHO* AREA* Vt lp**2) ) 

POINT THE DISTANCE BETWEEN THE ROTOR AXIS AND THE 

POINT WHERE THE AXIAL VELOCITIES ARE BEINC 
CALCULATED (FT) 

RAD THE ABSOLUTE VALUE OF POINT (FT) 

X THE RATIO OF VORTEX RINC DOWNSTREAM 

POSITION TO RINC RADIUS <Z( I )/RR( I ) > 

R THE RATIO OF POINT TO THE RING. RADIUS 

(RAD/RR ( I ) ) 


S.T.D.E.B.C, CIRC. .COMPONENTS OF TT1E EQUATION USED TO 


DETERMINE THE AXIAL VELOCITY DUE TO EACH 
RINC (FROM ‘AERODYNAMICS OF PROPULSION* . 
rUCHEMANN 8 WEBER. 1953) 

P.D INTERMEDIATE STORACEi VARI ABLES 

V(M.N) THE COMPONENT OF THE AXIAL VELOCITY AT 

POINT fl DUE TO RINC N (FT/SEC) 

VTOT ( 1 ) THE TOTAL AXIAL VELOCITY AT EACH OF 25 

POINTS DUEjTO SUMMING THE COMPONENTS FROM 
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C EACH VORTEX RING (PT/SEC) 

C VBAA THE AVERACt I NFLOW VELOCITT BETWEEN TVO 

c successive poiitts at which the axial 

C VELOCITIES WERE CALCULATED CFT/SEC) 

C DVHLDt THE DOWNLOAD COMPONENT FOR THE. AREA BETWEEN 

C TWO SUCCESSIVE POINTS AT WHICH THE AXIAL 

C VELOCITIES' WERE CALCULATED <L0> 

C DVHLD2 THE SUM Or»THE INDIVIDUAL DOWNLOAD 

C COMPONENTS*. I.e.. THE TOTAL DOWNLOAD (LB) 

C PRCNTDL THE PERCENT OF ROTOR » THRUST THAT IS DOWNLOAD 

C 

c 

CttwitttwMitwtwttwwimwtmmtm tintm imwm mma 

C 

C 


REAL r.KT(72> 

INTECER Y.TNR 

DIMENSION Z( 144 ) . COORD < 144 ) , F(36 ) ,RR( 144 ) , VT0T(7 1 ) ,V(7l ,144) 

C 

c 

c 

c 

C DUE TO THE FACT THAT THREE DIFFERENT ROTOR HEICHT TO ROTOR* 

C DIAMETER RATIOS WERE RUN WITH THE MODEL SYSTEM. THREE H/D 

C RATIOS CAN BE ANALYZED WITH THIS PROGRAM. TWO ARE IN GROUND 

C EFFECT (ICE) CONDITIONS. H/D-0.3 8 1.0. AND ONE IS AN OUT 

C OF GROUND EFFECT (OCE) CONDITION. H/D* I NFI NITY. THE SELECTION 

C TOR THE CROUND EFFECT CONDITION TO BE' ANALYZED IS MADE HERE. 

C 

C 

C 

C 


16 

VR ITE( 6,16) 

FORMAT ( * WHAT H/D RATIO WILL. BE USED? * 

14 

WR!TE<6, 14) 

FORMAT ( ’ (0.3.1 .0. OR 

OCE 

HOVER)?’ ) 

17 

VRITE(6. 17) 

FORMAT C (ENTER 0 FOR 

OCE 

HOVER) ’ ) 

18 

READ ( 3 , 18)B0VERD 
FORMAT ( C 1 2 . 3 ) 




c 

c 


OrtOTXtxx xxntxxxxxxxxxxxxxxxxxxMuxriiiCTXxirCTTTXOCTXXKXXJrctyyxr i r«»ixx : 

C 

C 

C NEXT. BASED ON A ROTOR SPEED OF 1000 RPM . THE TIME SPACINC OF 

C THE VORTEX RINGS IS CALCULATED USING THE RJNC SPACINC FACTOR (K>. 

C 

C 

C tttWttOtt mattmwWWMWCTMMtMMMIMWMXWMXXMMXW gCXXlXeX : 

c 

c 

K* .41946 
B* . 6837 
DO 20 I>1.72 

F ( I )*I*H/2.0 
KT( I )«IC*F( I ) 

20 CONTINUE 

C 
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c 


c 

G 

C A FIFTH ORDER I NTDIPOLATI KG POLYHUHI ALj IS USED TO DETERMINE THE 

C DOVTTSTREAM POSITION OF EACH VORTEX RING AS A FUNCTION OF TIME. 

C IF ONE OF THE ICE HOVER CONDITIONS IS BEING ANALYZED THE CALCULATION 

C IS TERMINATED AT THE RING THAT COINCIDES WITH THE GROUND PLANE AND 

C CONTROL IS TRANSFERRED TO STATEMENT 26. EACH VORTEX RING'S RADIUS 

C IS ALSO CALCULATED* AS A FUNCTION OF DOWNSTREAM POSITION. 

C 

C 

C 

c 

23 DO 40 1*1 .72 

zu )«3.342i*rr<i)**3.^i5;3269*inr<n**4.-24.eie7*icT( i > 

♦ **3 . + 12. 8294*KT ( I ) **2 . +7 . 2430*KT ( I ) 

RR<! >».66274*SQRT< I .✓< . 3*<2 . -EXF< -1 .6476*Z< I > ) ) > > # 

COORD ( I )«Z< D-. 123 

IF< (HOVERD. BQ .0 . 5 > -AND . (Z< 1 ) .CE. I . 133) .OR. 

♦ (BOVERD.EO. 1 .e).AND.<Z(l) .CE. 2 . 27 ) )TDEN 
RUHR! NC 8 I 

CO TO 26 
ENDIF 

40 CONTINUE 

TNR*72. 

CO TO 24 
C 
C 

C 

c 

C IN THE EVENT THAT ONE OF THE ICE ROVER CONDITIONS IS BE I NG ANALYZED 

C THIS SECTION FORMS ■ THE MIRROR IMAGE OF THE PREDETERMINED VORTEX 

C RING SYSTEM. IF THE OCE HOVER CONDITION IS BEINC ANALYZED THIS 

C SECTION IS BYPASSED AND CONTROL IS TRANSFERRED TO STATEMENT 24. 

C 

C 

C 

C 

26 Y-NUMRINCM 
Z<Y> «Z< KUKRI NC) 

RR(Y)«RR(NUMRINC) 

COORD ( Y )■ COORD ( NUHR I NC ) 

C 

C 

DO 27 1*2, HUMMING 
LOPE* I ♦NUMR1 NC- 
LIP*Y-| 

Z( LOPE ) *Z(Y)-MZ(Y)-Z<LIP) ) 

RR< LOPE) *RR< LIP) 

COORD ( LOPE )*Z< LOPE)-. 1979 

27 CONTINUE 
C 

c 

C 

C 

C THE VALUE OF THE TOTAL NUMBER OF RINCS <TNR> FOR ICE BOVERS . THRUST 
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C 

c 

c 

c 

c* 

c 

c 

24 


COEFF1C1EHT ICT). RIHC CIRCULATIOR (A) 1 ARD THE POSITIOH OF^YHE FIRST 
POIRT. VITH RESPECT TO THE ROTOR AXIS* (POUTT). ARE ROW OTTERED. 


TRR*2 * ® HUntt I RC 
CT* .0075 


A«2.226 

poirt»-i .3 
c 
c 


c 

c 

C TWO DO LOOPS ARE ROW USED TO CALCULATE EACH VORTEX RIRC’S AXIAL 

C VELOCITY COHPORTERT AT EACH P01HT OF IHTEREST. THE IRRER LOOP 

C CALCULATES THE VELOCITIES FROM EACH OF THE RIRCS ARD THE OUTER 

C LOOP IRCREflERTS THE POIRT AT WHICH THE VELOCITIES ARE CALCULATED. 

C 
C 

C 

C 


DO 63 H*1 .23 
C 

c 

DO 70 R«l ,TRR 

iF<POiirr.LT.e.> rad*-poirt 

IPIPOIRT.CT.O. ) RAD*POIRT 
X- -COORD* H)/RR<RJ 
R«RAD/RR<R) 

S=<4.*R>y<X**2.-MR+l . F**2. ) 

T* A LOCI 4. /SORT < \ .-S> ) 

U*TM IT- 1 . )/4. >«< I . -S )♦ ( 9 . /64 . )*<T- (7 . /6 . ) >* I I . -5)**2. 

E* I . + ( .3*<T-.5)*( I .-S> > + (3./l6. )*<T-( I3./12. ) >*< I .-S)**2- 
B=» I ./S0RT<X**2 . +IR+I . > **2 . ) 

+ . ) )/ <X**2 . -MR- 1 . )**2. ) ) ) 

D«U+C*E 

P-B*0 

A-ABS(A) 

IF( (BOVERDIRE. 0.00). ARD. (R.CE.Y) ) THER 

A* -A 
ERDIF 

Cl RC* Ay *2 . *3 . 141 56*RR( R ) ) 

V<H.R)»C1RC*P 
70 CORTIRUE 

C 
C 

P0IRT*P0IRT+.l23 
63 CORTIRUE 

C 
C 

> >■ 

C 

C 

C THE VELOCITY COrtPORTRTS DUE TO EACH OF 'THE VORTEX RIRCS ARE ROW 

C SimHED TO ARIVE AT A TOTAL AXIAL VELOCITY AT EACH OF THE 23 POIRTS 

C OF 1 RTEREST . 

C 



c 


jMmumw wni iinm u roiinmummHttwu iMW 


C 

C 

DO 73 J*1 .23 
C 

c 

do se i = i .Tim 

VTOTI J ) * VTOT( J ) +V ! J , I ) 

Be continue 

c 

c 

73 continue 

c 

c 

C 

C 

C THE TOTAL AXIAL VELOCITY AT EACH POINT OF INTEREST IS NOW PRINTED 

C OUT ALONG WITH THE ASSOCIATED VALUES OF ROTOR IIEICHT TO ROTOR 

C DI AJtETER RATIO. THRUST C0EFF1C 1 ENT . VORTEX RING SPAC1NC FACTOR AND 

C THE TOTAL NUMBER OF VORTEX RINCS USED IN THE VELOCITY CALCULATIONS. 

C 
C 

C 

c 

WRITE(6.33)CT.HOV£RD 

33 FORMAT ( • CT* * , F6 . 3 . ' H^D* * . F4 . 2 ) 

IF <HOVERD.EQ.e> THEN 

WRITE! 6 34) 

34 FORMAT ( * <H/D*0.ee SICNI FIES AH OCE HOVER CONDITION)*) 

ENDIF 

WR I TE ( 6 , 83 ) K . TNR 

83 FORMATUX.’ FOR K* ’ . FB . 7 . * AND’ .IX, 12.* VORTEX R I NCS ’ ) 

WRITE! 6 ,90 ) 

90 FORMAT! IX.' THE TOTAL AXIAL VELOCITY AT EACH POINT ALONC A. LINE ) 

WRITE!*, 95) 

95 FORMAT! IX.' CO I NC I DENTA L VI TH THE PLANE OF TDE W I NC I S : ' ) 

C 

C 

DO 97 J * I ,25 

WRITE (6,96) J.VTOTIJ) 

96 FORMAT ( / , 2X . 3HVT0T < . 12.2B)- , Cl 2 . 3 . 60JT/SEC ) 

97 CONTINUE 
C 

C 

C 

c 

C DOWNLOAD COMPONENTS ARE NOW COMPUTED BASED ON THE AVERACE AXIAL 

C VELOCITY BETWEEN TWO POINTS OF CONSIDERATION. THE TOTAL DOWN- 

C LOAD IS THEN COMPUTED AS THE i SUM OF THE COMPONENTS. 

C 

c 

C 

c 

DO 699 h j , 12 

VBARs (VTOT! I >*VTOT< I ♦ I ) )/2.0 

DVNLDl* . 5* . 80237 B*VBAR**2 .0* . 04896* I .98 
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wnnjM«BVfru»^wnn-Di 
w owrriiruE 

c 

c 

c 

c 

C THE FINAL RESULT FOR THE MODEL ROTOR SYSTEM THEORETICAL DOWNLOAD IS 

C HOW OUTPUT. 

C 

C 

C 

C 

WR I TE ( 6 . 390 > DVNLD2 : 

390 FORMAT < / . * THE MODEL ROTOR DOWNLOAD IS «',CI2.3.’ LB.’) 

PRCNTDL* < DVJ1LB2/0 .3)9100 
WRITER 6 .393) 

395 FORMAT** BASED ON A THRUST COEfTI Cl ENT! OF .0073 ’) 

VR I TE ( 6 . 400 ) PRCNTDL 

400 FORMAT ( ' THIS REPRESENTS' . IX. F3. 3..’ X OF ROTOR THRUST') 

C 

c 

STOP 

END 



